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ABSTRACT 


This study determines the performance and cost of four 10 MWe 
advanced solar thermal electric power plants sited in various regions 
of the continental United States. Each region has different insolation 
characteristics which result in varying collector field areas , plant 
performance, capital costs and energy costs. 

In the context of advanced technology, the solar plants are 
conceptualized to begin commercial operation in the year 2000. It is 
assumed that major subsystem performance will have improved 
substantially as compared to that of pilot plants currently operating 
or under construction. The net average annual system efficiency is 
therefore roughly twice that of current solar thermal electric power 
plant designs. Similarly, capital costs reflecting goals based on 
high-volume mass production that are considered to be appropriate for 
the year 2000 have been used. These costs, which are approximately an 
order of magnitude below the costs of current experimental projects, 
are believed to be achievable as a result of the anticipated sizeable 
solar penetration into the energy market in the 1990-2000 timeframe. 

The paraboloidal dish, central receiver, cylindrical parabolic 
trough, and compound parabolic concentrators comprise the advanced 
collector concepts studied. All concepts exhibit their best 
performance when sited in regional areas such as the sunbelt where the 
annual insolation is high. The regional variation in solar plant 
performance has been assessed in relation to the expected rise in the 
future cost of residential and commercial electricity supplied by 
conventional utility power systems in the same regions. 

The report contains a discussion of the regional insolation data 
base, a description of the solar systems performance and costs, and a 
presentation of a range for the forecast cost of conventional 
electricity by region and nationally over the next several decades. 


FOREWORD 


The advanced thermal technology work reported herein is a part of 
the thermal power systems activities of the Department of Energy’s 
Division of Solar Technology. A primary objective of this effort is to 
support development of advanced, low-cost, long-life, reliable solar 
thermal power systems which will supplement and eventually replace 
current fossil-fueled, electricity-generating plants. 

The National Aeronautics and Space Administration's (NASA) Jet 
Propulsion Laboratory (JPL) and Lewis Research Center (LeRC) were 
selected in 1977 to assist in managing and coordinating this work. 

These two organizations, working with universities, government 
agencies, industry and the scientific community in general, are to 
lead in developing new concepts and establishing a broad technology 
base in advanced small power systems which can be used to accelerate 
the commercialization of these systems. 

This report presents results of a study which determines the 
performance and cost characteristics of four selected advanced solar 
thermal power systems operating under widely varying insolation 
patterns associated with different regions of the country. The cost 
of electricity derived from the solar plants is compared to that from 
conventional (fossil and nuclear) power systems operating in the same 
region. 


The body of this report presents a summary of the key findings 
and conclusions regarding the potential attractiveness of advanced 
solar thermal power systems on a regional basis. Supporting detailed 
data and analyses are presented in a series of appendices covering 
(A) the regional insolation data base, (B) performance and cost 
characterization of the four selected solar thermal systems, including 
presentation in more detail than is provided in the body of the 
report, of results obtained by computer simulations of the power 
plants, and (C) a limited analysis of power plant performance and cost 
sensitivities to variations in certain major plant components and 
component costs. 
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STATEMENT OF THE PROBLEM AND ASSESSMENT OF THE RESULTS 


This study addresses the following two key questions: 

• What is the effect of regional insolation (solar 
energy) characteristics on the techno-economics of a 
spectrum of advanced solar thermal power plants? 

• How do the projected energy costs of solar thermal 
plants operating in different regions compare with 
energy costs of conventional systems in these 
regions? 

Since advanced solar thermal s stems in the 1990-2000 timeframe 
are assumed to achieve projected performance and cost goals, it is 
this timeframe that is considered here. The study entailed the 
selection and characterization of a set of advanced solar thermal 
power plants that are representative of the diverse spectrum of 
concepts. These selected power plant concepts were then configured to 
produce minimum energy costs at chosen sites within the continental 
United States. In this way, the effect of regional insolation 
characteristics on the design, performance, and economics of the 
selected set of solar thermal plants has been determined. Projections 
of energy costs for conventional systems are used as a basis for 
assessing the relative attractiveness of solar systems operating in 
different regions. 


A. SUMMARY 

A first-order screening to determine the feasibility of using 
solar thermal power as an energy source in different regions of the 
country has been accomplished, and the results of this screening are 
reported here. Whenever insolation was sufficient or whenever 
sufficient energy had been stored during prior periods of insolation 
to allow plant operation above some fraction of rated output, the 
solar plant was assumed to deliver power at partial to full rating. 

No plant was allowed to deliver power in excess of its rating at any 
time. The most significant conclusion of this study is that two-axis 

tracking solar thermal systems specifically, the central receiver 

and paraboloidal dish were relatively more efficient at all 

locations and thus, potentially attractive for electric power 
generation in all regions of the continental United States. 


B. REGIONAL INSOLATION CHARACTERISTICS 

The insolation data base for the continental United States, 
encompassing direct measurements as well as inferred data from 
meteorological observations, was analyzed to develop insolation 
contours. Solar insolation and the development of various models for 
this insolation are discussed in Appendix A of this report. These 
contours were then used as a guide in selecting eight sites 
representative of the wide range of regional insolation 
characteristics encountered within the continental United States. The 
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•elected sites end contours of annual mean or average daily direct 
normal insolation are shown in Figure 1. 

The chosen sites are compared in Table 1 in terms of four 
characteristic parameters which provide insight regarding some of the 
major differences among the sites. The locations and the year in which 
the insolation data were taken are given in descending order of annual 
total direct insolation. For concentrating solar collector systems, 
annual total direct insolation is a measure of the energy available. 

A large variation from a high of 2848 kWh/m 2 per year for Bars tow, 
California, to a low of 1253 kWh/m 2 per year for Maynard, 

Massachusetts, is shown. 

The year associated with the insolation data at each site was 
•elected by analyzing long-term insolation data at that site. A 
recent year, representative of long-term trends, was chosen since 
recent data are considered to be more reliable from the viewpoint of 
instrument accuracy and measuring/data reduction techniques. 

Referring to Table 1, it is seen that the annual total hours of 
insolation vary between 4000 h and 4500 h for all the sites except 
Maynard, MA, which has a value of only 2813 h. The annual average 
daily direct insolation is the annual total direct insolation divided 
by the number of days in the year. This parameter, which is 
proportional to total annual direct insolation, is used to identify 
contour curves in Figure 1. It is noted that the data in Table 1 and 
the contour curves of Figure 1 are in approximate agreement and that 
differences are indicative of uncertainties associated with developing 
contours based on limited data of varying quality. Microclimatic 


Table 1. Insolation Characteristics of Chosen Sites 


Local ion 
( year ) * 

Annual Total 
l) i rei l 1 n so 1 at i on 
kWh/(m J y) 

Annual Total 
Hours of 

Direct Insolation 
h/y 

Annual Average 
Daily Direct 
lnso lat ion 
kWh/lm^d) 

Annual Average 
Hourly Insolation 
(Based on Hours 
of Direct 
Insolation) 
kU/m 2 

Bars tow, OA 
(1976) 

784 8 

4494 

7.781 

0.6337 

Al buque r qut' , NM 
(197^) 

7* )8 

4 300 

6-9*) 

0. *907 

Omaha , N B 

1688 

4 76 3 

4.67* 

0.3960 

Medford, 0k 
(1973) 

1607 

4 376 

4 . 389 

' 0.3704 

Miami , Kl. 
(197*) 

14 74 

4138 

3.901 

0. 3447 

Cape Hatteras, NC 

(1973) 

1471 

4094 

3.893 

0. 3471 

Madison, U1 
(1971) 

1798 

4)43 

3.5*6 

0. 3132 

Maynard, MA 

1 7 Vi 

7813 

3.433 

0.44*6 


Year in which msol.it ion data was taken. 
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Figure 1. Insolation Patterns and Selected Sit 









effects associated with any chosen site add to the level of 
uncertainty and tend to require a large number of strategically 
located measuring stations to develop accurate contours. 

Nevertheless, the contour curves are valuable in providing an overall 
view of prevailing insolation patterns. 

In addition to the amount of available insolation, the intensity 
is also important with regard to the design and performance of solar 
thermal plants. Intensity level determines the efficiency of the 
receiver and engine when operating temperature is specified. In 
general, a higher intensity results in higher efficiency: For the 

receiver, re-radiation losses, which are essentially constant for a 
fixed temperature, become a smaller fraction of the heat processed by 
the receiver. As a result, receiver efficiency improves. If the 
engine rating is selected to match peak efficiency levels, the engine 
tends to operate closer to its design efficiency at near-design levels 
of solar insolation. A measure of insolation level is provided by the 
annual average hourly insolation as shown in Table i. In general, 
intensity level decreases as the annual total insolation or amount of 
insolation decreases. However, there are exceptions as exemplified by 
Maynard, MA, which has the lowest amount of insolation while having 
the third highest intensity- When insolation is available in Maynard 
MA, it tends to have high intensity whereas, in areas such as Miami, 
FL. atmospheric conditions dictate lower intensities over longer 
periods of time. 


C. SELECTED ADVANCED SOLAR THERMAL POWER PLANTS 

There is a spectrum of solar thermal power systems employing 
collectors ranging from low-temperature, fixed orientation to high- 
temperature, two-axis tracking concepts (Refs. 1 through 5). The four 
collector concepts illustrated in Figure 2 were chosen as being 
representative of the complete spectrum. 

o Compound Parabolic Concentrator (CPC) — a fixed- 

or ient at ion (non- tracking) collector which uses shaped 
reflecting surfaces to concentrate solar flux on a linear 
receiver; reflecting surface alignment is adjusted 
periodically (no more than twelve times oer year) for 
maximum performance. 

o Parabolic Trough — a one-axis tracking collector which 

utilizes a parabolic trough to reflect solar flux onto a 
linear receiver located along the trough's focal line; 
tracking is accomplished by rotating the surface about its 
focal line. 

o He liostat /Central Receive r — a two-axis tracking design 

which employs a field of reflecting surfaces, controlled 
so that the solar flux is beamed toward a tower-mounted 
receiver; a near-term embodiment of this concept is the 10 
MW pilot plant being constructed near Barstow, CA. 
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• Paraboloidal Dish — a modular two-axis tracking approach 
comprising a field of paraboloidal dish concentrators , 
each of which focuses solar flux into the aperture of a 
cavity receiver located at its focal point; a small heat 
engine-generator assembly is located behind the receiver 
at the focal point. 

A complete solar thermal power system was configured for each of 
the selected generic collector concepts using advanced technology 
projections for subsystems/components. These projections are 
generally reflective of the performance and cost goals for advanced 
solar systems that are supported by recent studies encompassing 
advanced design approaches and mass production manufacturing. These 
design and manufacturing studies indicate that the goals are within 
reach. 


The key factors governing the achievement of advanced solar 
thermal system energy cost targets (~50 mills/kWeh) are (1) the 
efficiency of the system in capturing and converting solar energy into 
electrical energy and (2) the unit cost of the collectors. These 
factors are interrelated. For a plant with a given output, the size 
or area of the collector field is inversely proportional to system 
efficiency. If the efficiency is low, the field becomes large and 
plant costs increase. To compensate, low-efficiency systems must 
employ low-cost collectors if they are to be competitive. 

In Table 2, the design performance characteristics of the four 
selected systems are presented. The projected performance levels are 
predicated on advancements in technology and refined designs for all 
key subsystems consistent with our projections for the existing 
conditions in the 1990-2000 timeframe. For the concentrator, 


Table 2. Design Performance Characteristics 

• Year 20^0 1 lartup • 10 MWe Plants 


Plant Type 

b* - h 

Central Receiver 

Trough 

CPC 

Operating Temperature, °C/ c F 

900/1650 

980/1800 

425/800 

225/440 

Concentrator 

0.89 

0.74 

0.74 


Tower Shading 

(1) 

0.98 

(1) 

(1) 

Atmospheric Attenuation 

(1) 

0.98 

(1) 

CD 

Receiver 

0.90 

0.68 

0.93 

_„<3) 

Concent r a tor /Receiver Efficiency 

(0.80) 

(0.63) 

(0.69) 

(0.41) 

Thermal Energy Transport 

0.99 (2) 

0.99 (2) 

0.92 

0.92 

Electrical Energy Transport 

0.96 

(1) 

(1) 

(1) 

Power Conversion 

0.47 

0.50 

0.36 

0.22 

Parasitic Power Losses 

0.93 

0.93 

0.93 

0.93 

NET P0WERPLANT SYSTEM EFFICIENCY 

0.34 

0.29 

0.21 

0.08 

(1) Not applicable. 

(2) Liquid metal, short length. 

(3) This data is calculated within 

the ANL coding. 
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back-silvered surfaces achieving reflectivity values of 0.90 to 0.95 
are employed. These estimates of reflectivity are based on the 
development of a thin, low-iron-content glass reflector technology, 
such as the fused glass process. Manufacturing techniques providing 
surfaces with small slope errors 50.1° are assumed. Receiver 
losses are minimized by employing high-quality insulating materials. 

Each of tie different systems has an optimum operating 
temperature rang 2. This range is determined via a tradeoff between 
receiver efficiency, which decreases with temperature (due primarily 
to re-radiation losses), and engine efficiency, which increases with 
temperature (Ref. 3). Two-axis tracking systems are able to obtain 
higher concentration ratios and correspondingly smaller receiver 
areas. In turn, these smaller areas permit the achievement of higher 
temperatures before re-radiation losses become dominant. Thus, as 
seen from Table 2, two-axis tracking systems have the highest design 
operating temperatures (in the range of 1500°F to 2000°F for the 
selected paraboloidal dish and central receiver systems). The 
one-axis tracking trough and fixed-orientation CPC attain lower 
concentration ratios and the selected designs for these configurations 
operate at lower temperatures of 800°F and 440°F, respectively. 

The temperature level determines the theoretical upper-bound 
Carnot efficiency for the energy conversion system. Present heat 
engines achieve approximately 60% of Carnot, whereas the very advanced 
systems are projected to achieve values of 80% Carnot (Ref. 3). For 
the 1990-2000 timeframe, an advanced technology target of 70% Carnot 
was deemed to be appropriate, and this value was used for all systems. 

Stirling engines were assumed for the paraboloidal dish system. 
Near-term development of these engines is focused on designs operating 
at ~ 1500°F. A modest increase to 1650°F is projected for the 
1990-2000 timeframe. An advanced Bray ton/Rankine combined cycle was 
projected for the central receiver. To achieve the efficiency 
benefits of this particular combined cycle, high temperatures of 
~1800°F are desirable and 1800°F was used in the study. The 
parabolic trough employs either steam or organic Rankine power 
conversion. Advanced organic Rankine engines are best suited to the 
CPC. 


Development of Stirling engines is concentrated on small sizes 
suitable for the paraboloidal dish system. Small Brayton engines are 
also candidates for dish systems since earlier studies (Ref. 2) 
indicate that performance reasonably close to that of the Stirling 
system can be achieved. For the large ( ~ lOMWe) engines needed for 
the central receiver concept, the Brayton/Rankine system was selected 
since these types of engines are available in this size range. Large 
Stirling engines, somewhat similar to large multi-piston diesel 
engines, are also candidates, but as mentioned above, current 
developmental activity is focused on the smaller sizes. 


ft 
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It is noted from Table 2 that power conversion efficiency 
decreases as operating temperature decreases. This follows the 
expected trend inferred from the Carnot efficiency variation with 
temperature. The power conversion efficiency also includes the effect 
of mechanical/gearbox and electric generator efficiencies. Lower 
mechanical-generator efficiencies are assumed for the small power 
units that are a part of the modular paraboloidal dish concept than 
are assumed for the other chosen concepts, which employ large 
centralized energy conversion units. 

In terms of net or overall system efficiency, the optically more 
efficient two-axis tracking concepts achieve higher values than one- 
axis and fixed-orientation systems. However, it is emphasized that 
this greater efficiency is achieved at the price of greater complexity 
in the articulation of concentrating or reflecting surfaces. The 
paraboloidal dish appears to have the highest potential for 
efficiency, primarily due to its optical/tracking characteristics 
which avoid the cosine losses encountered to differing degrees by the 
other systems. Basically, the dish tracks so that the collector 
aperture is normal to the solar beam radiation whereas, for the other 
systems, the collector aperture cannot be maintained normal to the 
beam over the entire period of insolation availability. 

Fixed-orientation systems, as represented by the CPC, employ 
relatively lew concentration ratios of the order of five. The 
receiver is located at the base of the reflecting surfaces, as seen in 
Figu'.e 2, and therefore can also absorb a portion of the diffuse 
radiation. The other systems only accept the direct beam radiation. 
When the direct beam radiation is blocked by cloud cover, these other 
systems cannot function. Systems, such as the CPC, which use some 
diffuse energy may still be able to function but, generally, when beam 
radiation is blocked, the entire level of insolation is low, and 
performance is correspondingly low. 

Systems which convert diffuse solar energy could be advantageous 
for applications involving low-temperature heat, particularly in 
cloudy regions. When considering only power generation as in the 
present study, this advantage is greatly diminished. Although systems 
such as the CPC appear low in terms of performance potential in the 
context of the present study, it is emphasized that there is a wide 
range of solar thermal applications and that the CPC could be a strong 
candidate for some of these applications. Additionally, the CPC is at 
a more advanced stage of development and can be readily implemented 
for near-term applications. 

As an additional point, it is noted in Table 2 that the net 
system efficiency applies to operation at design conditions where none 
of the energy passes through storage. To operate at high capacity 
factors (i.e., to deliver rated power for periods beyond the time of 
insolation availability), it is necessary to introduce energy 
storage. For the dish-Stirling system, advanced battery storage 
having a throughput efficiency of 75% is assumed. The other systems 
assume thermal storage with a throughput efficiency of 85%. The 
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effect of storage efficiency on total system performance is clearly a 
function of the amount of energy that passes through storage* The 
amount of energy passing through storage is determined by the desired 
capacity factor and the characteristics of the insolation. Thus, 
regional variations in insolation patterns will affect storage 
requirements and system performance* 

The direct cost data for each of the selected systems is 
presented in Table 3. As determined from sensitivity studies, the key 
cost driver is unit concentrator cost. Since solar energy has 
relatively low intensity at the earth's surface, it is necessary to 
use collector fields of large area to capture and concentrate 
sufficient energy for applications such as electrical power generation. 
Since large areas are required, the cost per unit area must be small in 
order to achieve cost-competitive systems. 

Projections of low-cosf, advanced concentrators, as shown in 
Table 3, are predicated on high-volume mass production (~ It)** 
units/year) of advanced designs specifically tailored for mass 
production. Most of the recent studies (e.g., Refs 6 and 7) concerned 
with high-volume mass production have treated the concentrator or 
heliostat of the central receiver system. Studies have encompassed 
the heliostat designed for the 10 MWe pilot plant being constructed 
near Barstow, CA, as well as advanced second-generation concepts. 

These studies indicate that a unit cost of $65/m 2 (as shown in Table 
3) is achievable for a completely installed system. 


Table 3. Solar Subsystems Direct (D Cost Data 


• 10 MWe Power Plants 

• Vear 2000 Startup 

• 1979 Base Year 

Dollars 


Major Subsystems 


Paraboloidal 

Dish 

Central 
Rece i ver 

Parabo 1 ic 
Trough 

Compound 

Parabolic 

Concentrator 

Operating Temperature c C/°F 


900/1650 

980/1800 

425/800 

225/440 

Land, $/acre ($1000 - $5000) /acre 


__(2) 

_„(2) 

(2) 

m 

Concentrators, $/i n*-’ 


83 

65 

) 

) 





)70< 3 > 

)50< 4 > 

Cavity Receivers, S/m^^ 


10 

) 

) 

> 




) 1 75 $/kWeW 



Energy Transport 


40 $/kWe 

) 

60 $/kWt 

50 $/kWt 

Storage St rue lures/Med i um^ 


30 $/kWe-h 

26 $/kW-h 

24 $/kW-h 

21 $/kW- h 

Heat Engine Equipment, $/kWe 


60 

60 

60 

60 

Balance of Plant, $/kWe 


85 

85 

85 

oo 

Operations and Maintenance Costs, 

7«0 

1 .0 

1.0 

1 .0 

1.0 


(1) Direct costs do not. include: spares/curit ingeiicy (52), indirect costs (10%), or interest /during 

construct ion .* 

(2 ) Land costs are ite dependent. 

(3), (4) Estimated mass production costs: JPL, coordinated with Sandia, Albuquerque and Argonne National 
Laboratory . 

(5) Receiver cost is normalized to concent r.i tor aperture area. 

(6) Includes tower, structure, receiver, and short piping transport in tower mounted power conversion 
sy st cm modu 1 e . 

(7) Based on 5-6 h ot storage. 

(8) Operations and maintenance (O&M) costs for the first year are assumed to be a percentage of 
the total capital investment cost;06M costs escalate at the rate of IT per vear. 
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The values for unit cost shown in Table 3 were determined by 
using these studies as a point of reference in conjunction with goals 
established for the other systems. Generally, the paraboloidal dish 
employs doubly curved surfaces and is structurally more complex than 
the heliostat. Hence, it was assigned a unit cost of $85/m*, which 
is higher than the heliostat and within the established cost target 
range of 70-100 $/m 2 . Additionally, recent JPL contracted studies 
for dish systems projected installed costs in the 70-120 $/m^ range 
based on a mass production rate of 100,000 units per year. 

Systems such as the one-axis tracking trough and the CPC were 
judged to be simpler and less costly on a unit area basis. 

Laboratories directing the development of these systems, Sandia 
Laboratories, Albuquerque, for the trough, and Argonne National 
Laboratory for the CPC, concurred with the judgment that troughs and 
CPC’s should have lower costs than two-axis tracking heliostats and 
dishes. Since high-volume mass production studies comparable to those 
made for the heliostat have not been performed, there is uncertainty 
regarding the absolute magnitude but there is general agreement 
regarding the trends. 

Storage costs shown in Table 3 correspond to advanced technology 
projections. Since Department of Energy cost targets are in the 
process of being updated, sensitivities to the assumed values are 
presented in Appendix C. Costs for the heat engine and associated 
conversion equipment are values based on large-volume mass production 
of small engines via techniques developed within the automotive 
industry. Thus, the costs shown in Table 3 of $60/kWe strictly 
pertain to projections for the small heat engine employed in the 
paraboloidal dish concept. However, if the other systems are used for 
small power applications and if sizeable penetrations into the power 
market are made, some level of mass production manufacturing could 
also occur in their cases. Thus, it was decided to use the same 
energy conversion system unit cost for all systems, because it is 
believed that the equivalent unit costs of large engines could, at 
best, approach the mass production unit costs of the small engines. 

Balance-of-plant and operations-and-maintenance costs (O&M) were 
based on detailed studies for the central receiver and estimates by 
the Electric Power Research Institute (Ref 8). These factors depend 
on detailed design considerations. Since detailed designs do not 
exist for the projected advanced systems, only approximate estimates 
could be made. Therefore, it was decided to use the same factors for 
all the systems. 

As additional support for the use of the same factors, it is 
noted that there is no clear way to distinguish among the different 
systems without detailed study. For example, two-axis tracking 
systems are generally more complex than one-axis and fixed-orientation 

systems and would probably require more maintenance per unit area of 

collector field. However, the simpler but less-efficient systems will 

require a larger collector area to deliver the same power. These 

effects tend to be off-setting. 
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D. 


EFFECT OF REGIONAL INSOLATION ON SOLAR POWER PLANTS 


Regional insolation characteristics have a significant impact on 
the design, performance and cost of solar thermal power plants. For a 
given plant output, the amount and intensity of the available 
insolation at a site are major factors in determining the size of the 
collector field and, hence, its cost. As noted previously, the 
collector field is the key cost driver for solar plants. 

Based on performance potential, the southwest sunbelt is the 
best location for solar plants because it provides the greatest amount 
and intensity of insolation. However, depending on the performance 
level that they can achieve with prevailing insolation 
characteristics, there may be many potential applications for solar 
thermal systems outside the sunbelt. 

The effect of insolation on solar power plant design and techno- 
economics was determined by using a previously developed performance/ 
economics computer program (Ref. 9). The primary input of this 
program is hourly insolation data as supplied by SOLMET tapes. The 
performance portion of the program takes all the parametric 
characteristics of the power plant subsystems and links them into a 
simulation of a total system which operates according to specified 
control strategies; these strategies involve delivery of energy, 
storage of energy, and discharge of storage. Basically, the 
performance simulation determines the fraction of available insolation 
that is delivered as electrical energy. 

Also, the required collector field area and energy storage 
requirements to deliver a specified level of energy are determined 
from the performance simulation. Component/subsystem unit costs as 
given in Table 3 are used to determine corresponding capital costs. 

The well-established economic methodology for determining levelized 
bus-bar energy costs for utility-owned power plants is then employed 
(Ref. B-10). A computer optimization technique is employed to 
delineate power plant configurations which provide minimum energy 
costs for a specified level of energy delivery (i.e. capacity 
factor). Further details are provided in Appendix B. 

A set of three charts (Figures 3a, b, and c) summarizes the 
effect of regional insolation characteristics on solar power plant 
energy costs. In addition to expected overall trends of higher energy 
costs in regions of low insolation, some of the more subtle influences 
of differing regional insolation characteristics are delineated. The 
cnarts show energy cost as a function of annual average daily direct 
insolation (which is proportional to total annual direct insolation) 
for the selected advanced solar thermal power plants. The set of 
three charts is used to illustrate the effect of capacity factor 
(i.e., Figure 3a is for the capacity factor corresponding to zero 
storage whereas Figures 3b and 3c pertain to capacity factors of 0.4 
and 0.53, respectively). For the latter two higher capacity factors, 
storage is required. The effect of regional insolation 
characteristics on storage and its influence on plant techno-economics 
is discerned by comparing the three charts. 
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Figure 3a. Performance of Three Types of Solar Thermal Power 
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Figure 3b. Performance of Three Types of Solar Thermal Power 
Plants at Various Geographic Locations 
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When storage is not used, the collector field is sized to 
produce rated power during peak insolation periods and there is 
essentially no energy that is wasted. As shown in Figure 3a, the 
energy costs are lowest in the sunbelt region as represented by 
Albuquerque, New Mexico, and Barstow, California. As the annual 
average daily direct insolation decreases, the levelized cost of 
electricity or energy cost increases. 

There are particular perturbations about the overall trend that 
are caused by site-specific insolation characteristics. First, it is 
noted that, for the selected nominal designs, the central receiver has 
a slightly higher energy cost than the paraboloidal-dish-Stirling in 
the sunbelt and slightly lower costs in locations such as Medford, 
Oregon, and Maynard, Massachusetts, where the insolation level is 
lower. This cross-over is explained in terms of differing off-load 
power conversion efficiency characteristics for the central receiver 
with combined Brayton/Rankine power conversion as compared to the 
parabolic dish with an advanced Stirling heat engine. 

Closed-cycle Brayton engines typically have off-load curves with 
peak efficiency at 50-60X of rated load. This type of characteristic 
is therefore evident in the combined cycle system. The advanced 
Stirling engine used in the nominal designs achieves peak efficiency 
at rated power. Further, the engine located at the focal point of 
each dish was sized to handle peak insolation levels of the order of 
I kW/m? at all sites. 

When a sizeable fraction of the available energy is furnished at 
near the peak insolation value, as in the sunbelt, engines with 
off-load curves having highest efficiencies near the rated peak design 
point provide higher performance. However, where a greater portion of 
the energy occurs at lower intensity levels, as in Maynard, MA, or 
Medford, OR, off-load characteristics with maximum efficiencies 
occurring at partial load tend to be favored. This is one reason why 
the central receiver with Brayton/Rankine improves relative to the 
dish-Stirling in regions of low insolation. 

The perturbations for the parabolic trough (as seen in Figure 
3a) are also related to site-specific insolation characteristics which 
affect system performance by causing the engine to operate at 
different off-load settings for different amounts of time. The 
parabolic trough employs a Rankine conversion jystem with off-load 
characteristics similar to the Stirling engine (i.e., maximum 
efficiency at rated power). 

The solar plant design can be perturbed and re-optimized to 
accomodate the perturbations and site-specific insolation 
characteristics noted above. For example, engine off-load curves 
could be altered through changes in engine design, (i.e., engine 
designs could be tailored to meet desired off-load characteristics). 
The design rating of the engine could be selected to match 
lower-than-peak insolation conditions in regions where the average 
insolation is low, and excess energy could be wasted or stored during 
the few peak insolation periods. 
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An attempt at such re-optimization was not made in the present 
study because re-optimization to tailor plants to site-specific regional 
insolation characteristics involves complex and detailed design trade- 
offs, Therefore y as a caveat, it is noted that values of energy cost 
lover than those shown could be achieved by further design optimizations. 

Comparisons of Figures 3b and 3c with 3a reveal the influence of 
storage. First, the overall magnitude of energy costs shifts to 
slightly higher values as capacity factor and, hence, storage increase. 
This shift is attributed primarily to added capital costs for storage 
as well as the loss in energy associated with passing through storage. 

For the chosen nominal designs, the paraboloidal dish system is penal- 
ized relative to the other systems since an advanced battery storage 
throughput efficiency of 75Z was used as compared to 85Z for the other 
systems employing advanced thermal storage. Thf re are indications that 
advanced battery systems such as the sodium-sulfur system may achieve 
efficiencies higher than 75Z and thereby reduce differences between 
battery and thermal storage. The values used in this study are nomi- 
nal projections, and results should be interpreted in this context. 

The difference between the paraboloidal dish and central receiver 
is slightly greater at a capacity factor of CF « 0.4 (Figure 3b) as 
compared to a capacity factor of CF * 0.55 (Figure 3c). Here, higher 
capital cos ts for thermal storage systems are partially offsetting 
their efficiency advantages over battery systems. The effect of uncer- 
tainties in nominal storage cost values are presented in Appendix C 
because recent activities at NASA/LeRC have indicated lower thermal 
storage cost goals than were used in the present study. 

In addition to an upward shift in energy cost levels with the 
addition of storage, it is seen that site-specific perturbations are 
greater. This trend is due to insolation characteristics which impact 
the storage system in different ways. Basically, the most desirable 
insolation pattern would be a cyclic (diurnal) variation that persists 
essentially unchanged over the entire year. The storage system could 
then be sized and optimized to meet this variation. However, seasonal 
variations in cyclic insolation characteristics can be very large, par- 
ticularly for some regions outside the sunbelt. For example, there is 
a strong seasonal variation in Medford, OR, and it is seen that pertur- 
bations for Medford become very large at a capacity factor of 0.55 
(see Figure 3c). 

To achieve higher capacity factors, the collector field area is 
increased and the excess energy is stored and delivered during even- 
ings or periods of insolation non-availability. Effective use of stor- 
age is difficult if the seasonal variations are such that there is one 
prolonged period of strong insolation and another prolonged period of 
poor insolation. During the period that strong insolation having reg- 
ular diurnal insolation (cycles) is available, the field and storage 
size can be increased to the point where rated power is delivered dur- 
ing the full evening period. However, since this type of operation 
can occur for only a part of the year, the increase in capacity factor 
(based on annual energy delivered) will be proportional to the portion 
of the year for which strong insolation occurs. Storage will be 
largely under-utilized during the prolonged period of low insolation 
availability. 
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Depending upon seasonal inaolation characteristics, there is a 
capacity factor corresponding to the condition where the energy during the 
strong insolation period is maximally utilised. If one attempts to 
achieve even higher capacity factors by further increases in collector 
field area and storage, there will be insufficient time to deliver all the 
stored energy by the next sunrise. To accomodate this buildup of excess 
energy, large seasonal storage must be supplied and/or energy must be 
wasted. Costs for long-duration seasonal storage become very large for 
the assumed thermal and battery storage systems and the computer 
optimization algorithm results in wasted energy. This results in a rapid 
increase in energy costs. Thus, sites such as Medford, OR, having large 
seasonal variations in insolation characteristics will generally have a 
lower practical operating range of capacity factors. 

These effects of seasonal variations in insolation are illustra- 
ted in Figure 4, where the relationship between collector area and 
capacity factor is shown for the paraboloidal dish system. Other systems 
show similar trends. For low capacity factors, the relation between 
collector area and capacity factor is essentially linear. This trend is 
reasonable since capacity factor is a measure of the energy delivered and 
the energy delivered is directly proportional to the energy collected and 
to the collector crea. 

Also, in the regions of low capacity factors, the sites with less 
total annual direct insolation (see Table 1) require more area to deliver 
a given quantity of energy. Here, a given quantity of energy delivery is 
denoted by a fixed value of capacity factor. This effect is apparent in 
the steep slopes for curves corresponding to sites having low insolation. 

Departures from the linear relationship occur at higher capacity 
factors when storage is introduced. Since a portion of the collected 
energy is lost in storage, requirements for collector area increase and 
the slopes of the curves tend to be steeper as the level of storage 
increases. The rate at which the curves increase in slope is strongly 
affected by seasonal variations. For example, the slope for Medford, OR, 
increases sharply after a capacity factor of 0.4. The large area and 
storage at a capacity factor of 0.55 explains the reason for the large 
difference in energy cost for Medford, OR, as compared with other sites 
shown in Figure 3c. 

It is noted that the fixed-orientation system, as represented by the 
CPC, is presented only on Figure 3c and only to show its general magnitude 
of energy costs relative to the other system. Generally it exhibits 
trends similar to the other systems. Its ability to use a portion of the 
diffuse energy does not provide a significant advantage at the operating 
temperature levels deemed necessary for power generation. If a given 
temperature level is maintained when insolation is low (large diffuse 
component), re-radiation losses from the receiver become a large fraction 
of the incoming energy and system efficiency drops. Engine performance 
suffers if a lower operating temperature is used to reduce re-radiation 
losses. Thus, it is difficult to generate power in any efficient way when 
insolation is low and in diffuse form. Diffuse solar energy may be more 
effectively utilized in low-temperature thermal applications. 

Appendix B contains supporting detailed analyses which corroborate 
the findings discussed above. 
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REGIONAL COMPARISON OF SOLAR THERMAL WITH CONVENTIONAL ENERGY 



E. 

At compared to performance in the solar-intensive Southwest sun- 
belt, the performance of solar plants will be degraded in regions such as 
the Northeast. However, the cost of electrical energy derived from 
conventional sources presently shows strong regional variations. 

Forecasts such as those of Data Resources Inc. (Ref 9), indicate that 
relatively strong regional differences will persist through at least the 
year 2000. 

Moreover, regarding electrical costs from coal and gas turbine plants, 
there are large uncertainties that arise from potentially large variations 
in fuel prices coupled with capital cost increases that depend on the 
stringency of environmental control policies. Nuclear power also faces an 
uncertain future with regard to the nature and extent of the steps 
necessary to gain public acceptance for its expanded usage. 

The projected energy costs for solar thermal plants with no storage 
are compared in Figure 5 with energy cost forecasts for conventional power 
plants. The ranges shown for gas turbine and coal plants were taken from 
the work of Anderson, Bowers, et al.(Ref. 12) of Oak Ridge National 
Laboratory. The basis of their projections is as follows: 

Fuel Cost Range Capacity Factor 

Gas Turbine Residual Oil 14 to 21 $/bbl 0.3 

Coal Plant Low Sulfur Coal 28 to 49 $/ton 0.6 

Independent studies by JPL (Ref. 13) closelv agree with the range 
determined by Anderson and Bowers. 

The Data Resources Inc. (DRl) forecast (Ref. 11) is also shown on 
Figure S. The band corresponds to worst and best case scenarios 
formulated by DRI. Since nuclear power was projected to have the lowest 
costs, the DRI forecast is predicated on a strong shift to nuclear power 
based on an econometric supply-demand model. As noted previously, the 
occurrence of a strong shift to nuclear power is presently a highly 
uncertain proposition. 

For solar plants with no storage, which can be considered as the 
solar counterpart of the peaking gas turbine plants, energy costs equal to 
or less than the DRI forecast are projected for the central receiver and 
paraboloidal dish systems in the sunbelt region. For regions outside of 
the sunbelt, solar plant energy costs f-'i the central receiver and dish 
are within about 40% of the DRI foreca-'. It is noted that the DRI 
forecast is for a mix of plants including old and new plants. New plants 
would tend to be more costly and have higher energy costs. When 
considering the implementation of solar plants, their energy costs should 
be compared with the energy costs of new conventional plants to be 
constructed in the same time period. Thus, a strictly one-to-one 
comparison with the DRI forecast cannot be made. Rather, the DRI forecast 
is regarded as a general reference level which is of pa icular value in 
illustrating regional variations in energy costs for conventional 
svstems. Also, it is noted that no-storage solar plants are particularly 
suited for summer-peaking electrical demands where peak insolation ind 
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Figure 4. Collector Area as a Function of Capacity Factor 
for the 10 MWe Paraboloidal Dish Stirling 
Solar Thermal Power Plant at Various Sites 
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peak demand are coincident* The matching of insolation 
characteristics with load demand curves and the integration of the 
solar plant into utility grid systems were beyond the scope of this 
study. Instead, this study focused on determining energy costs for a 
simple fixed-demand characteristic. Energy costs determined in this 
manner are regarded as a measure of potential attractiveness. 

From Figure 5, it can be inferred that paraboloidal dish and 
central receiver solar thermal peaking (no storage) plants will 
achieve energy costs that are competitive with residual-oil, 
gas-turbine peaking plants for virtually all sites within the 
continental United States. The parabolic trough will be competitive 
in the sunbelt. To reach this competitive status will require 
advanced technology development, low-cost designs suited to mass 
production, and sufficient market penetration to achieve the 
high-volume production associated with low costs. 


F. CONCLUSIONS AND RECOMMENDATIONS 

The primary conclusion from this study is that the two-axis 
tracking solar thermal systems (namely the central receiver and 
paraboloidal dish systems) are potentially attractive for electric 
power generation in virtually all regions of the continental United 
States and the relative advantages of these systems over the other 
configurations hold for all regions. The most promising region is the 
Southwest sunbelt. However, in the context of large uncertainties 
regarding future energy costs from fossil and nuclear power plants, 
solar thermal power generation is potentially viable, even in the 
Northeast . 

The parabolic trough could be competitive in the sunbelt for 
power generation. The fixed-orientation system does not appear to be 
competitive for strictly power generation applications. However, it 
is strongly emphasized that solar thermal systems are suitable for a 
wide range of applications requiring heat at different temperature 
levels. Thus, systems such as the trough could be used to supply both 
heat and electricity for total energy applications (the concept being 
pursued at Sandia Laboratories, Albuquerque), and the CPC could be 
developed to supply low-temperature heat. The present study does not 
address these possibilities, and the results should be interpreted in 
this context. 

The present study is regarded as a first-order screening to 
assess the relative viability of using solar thermal power in 
different regions of the nation. For this purpose, a simple constant- 
demand load was assumed, where the solar plant delivered power during 
periods of insolation availability at levels equal to or less than 
rated power. When solar energy levels exceeded that needed to supply 
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plant rating, energy was stored and delivered at rated power during 
periods when insolation was not available. If excess energy was 
available and storage was full, energy was wasted. It is recognized 
that there are strong site-dependent differences in load 
characteristics that will affect plant design, operation and energy 
cos ts * It is therefore recommend ed that following studies incorporate 
the effect of differences in regional demand characteristics and focus 
mainly on the two-axis tracking systems which have been shown from 
this study to be potentially attractive in all regions* 

It is also recommended that a wide spectrum of solar thermal 
applications involving supply of heat over a range of temperatures be 
examined in the context of variations in regional insolation 
patterns. This is particularly important for systems requiring energy 
storage since it was shown that strong seasonal variations in 
insolation patterns have a severe impact on storage system 
requirements which, in turn, exert a major influence on plant 
operation and economics. 
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I. 


INTRODUCTION 


The evaluation of different thermal power plant design con- 
cepts requires a knowledge of the solar insolation available across the 
continental U.S. This appendix sumnnrizes the various aspects of 
insolation that must be considered in such evaluations, analyzes the 
current status of the data relating to these aspects with respect to 
accuracy and completeness, and discusses the basic information 
necessary to compare the technical and economic feasibility of various 
solar thermal power systems. 


II. THEORY OF INSOLATION 
A. BASIC CONCEPTS 

In addition to the actual magnitude of the insolation values, 
the two aspects of insolation which are of major importance in 
evaluating solar thermal power systems are its directional and 
temporal properties. 

The radiance of the solar radiation field at a given point in 
space, R(tl), represents the power per unit area normal to a specified 
direction, per unit solid angle, and is a function of the direction 
specified. Outside the earth's atmosphere this radiance has a value of 

R(0) - 1.988 x 10 4 ~2 — 

ra steradian 

in the direction of the solar disc, and is essentially zero for other 
directions.* This value is independent of the sun-earth distance and 
does not therefore have an annual time dependence. Being a measure of 
intrinsic solar properties it is also believed to be constant in time 
over very long periods, although recent studies have suggested 
possible variations of a few percent over the past few decades. 

The solar constant is obtained from the above radiance by inte- 
grating over the solar disc at a sun-earth distance of one astronomical 
unit and has the value of 1.354 kW/m^. This value is also constant 
in time and does not depend upon direction. 

In contrast to the simple conditions outside the earth's atmos- 
phere, scattering and absorption phenomena result in a much more 
complex situation for locations at or near the earth's surface. 

These effects may be summarized by comparing the radiance at the 
earth's surface with that outside the earth's atmosphere, and such a 
comparison is shown in Fig. A-l . It is seen that the terrestrial plot 
of R(u) differs from the exoatmospher ic plot in having a smaller value 
across the solar disc ( 1.5 x 10^ kw/ra^ - steradian) and non-zero 
valufco for angles away from the disc. 


*This is an average value across the solar disc and ignores limb 
darkening. 
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The above discussion has dealt with the total power integrated 
over all wavelengths* Not shown in the comparison of Fig* A-l is the 
fact that scattering and absorption produce a wavelength dependence 
(The sun is yellow; the sky is blue). Also not shown is the fact that 
the plot has a strong temporal dependence as well as a dependence upon 
the location on the earth's surface where the measurements are made* 

The temporal and spatial variations of R(0) are of great importance to 
solar thermal power system operation and will be discussed in greater 
detail in Section A-III. 

B. INSOLATION MEASUREMENTS 

Radiation measuring instruments (radiometers ) are designed to 
measure power levels integrated over some fixed area and solid angle 
determined by the instrument geometry. The quantity they specify is 
the power per unit area normal to a specified direction, i.e., the 
irradiance over the given solid angle. 

Insolation instruments differ in the solid angle used and the 
orientation of the instrument axis relative to the sun's center. They 
fall into two categories depending upon the values of these parameters* 

The pyrheliometer has a relatively small field angle (typically 
5°-15°), and is maintained with its axis directed towards the sun 
by means of a clock-driven equatorial mount. The quantity thus 
measured is called the direct normal insolation and is related to the 
radiance shown in Fig* A-l by the following expression, 

I DN " [q R(0) dfi , 

where is the solid angle corresponding to the instrument's field 
angle. 

The pyranometer has a 180° field angle and hence measures the 
insolation over a solid angle of 2tt steradians. The orientation of 
the instrument axis is fixed relative to the local vertical direction 
and is usually parallel to this direction. In this case, the 
resulting measure is called the total horizontal insolation or total 
hemispheric insolation , I^h* Since the direction to the sun's 
center relative to the local vertical direction has a complex temporal 
dependence which, in turn, depends upon the geographical location, the 
relationship between I^h and R(6) is no longer simple. 

The pyranometer is sometimes used with a shade ring which is 
mounted ru as to prevent the direct radiation from the sun from 
reaching the instrument's entrance pupil. In this case the instrument 
measures the diffuse horizontal insolation , I^h* 

As can be seen from Fig. A-l, a pyrheliometer with a total field 
angle of several degrees will measure not only the direct radiation 
from the sun (field angle « 1/2°) , but also a certain amount of 
diffuse radiation coming from the sky. While generally small, this 
diffuse radiation which enters the pyrheliometer (called the circum- 
solar insolation, t cs ) can be significant in the presence of strong 
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atmospheric scattering by water vapor. The dashed curve in Fig. A-l 
corresponds to a hazy atmospheric condition and illustrates the strong 
atmospheric scattering power "f water vapor and submi cron-s ize water 
droplets. 

If the irradiance measured by a pyrheliometer with a 1/2° 
field angle is defined as the true direct normal component, IdnO > 
then the quantities defined above are related as follows: 

X DN " *DNO ♦ ^s* 

I TH " J DN co * z ♦ x dH 

“ ^DNO + I cs^ C0B 1 * x dH» 
where: Z * zenith angle of sun. 


III. TEMPORAL AND SPATIAL VARIATIONS OF THE INSOLATION OBSERVED AT 

THE EARTH* S SURFACE 

There are a number of different phenomena which produce temporal 
and spatial variations in insolation, and these will be treated 
individually in the following sections. 

A. EARTH'S ORBITAL ECCENTRIC ITT 

The earth-sun distance varies between 1.47 x 10® km (peri- 
helion - Jan. 3) and 1.52 x 10® km (aphelion - July 5) causing the 
apparent diameter of the solar disc to vary from 32.55' to 31.48*. A 
focusing collector with a field angle slightly larger than the former 
value, aimed at the disc's center will thus detect an annual variance 
in irradiance of about +3.4% from this cause. 

B. DIURNAL VARIATION 

In the absence of an atmosphere, a focusing collector continu- 
ously aimed at the solar disc wou^d measure a nearly rectangular time 
dependence for the irradiance from sunrise to sunset. The effect of 
the atmosphere is twofold. First, atmospheric refraction extends the 
length of a day by about 4.5 minutes in the continental U.S., and 
second, atmospheric absorption and scattering round off the 
rectangular shape to produce typical clear-day curves such as those 
shown in Fig. A-2 for the direct-normal component. 

C. LATITUDE AND SEASONAL EFFECTS 

In the absence of atmospheric effects, the direct normal i* rad- 
iance integrated over an entire year would be the same anyv ;re on 
earth, neglecting the aforementioned orbital distance effect. Lcn^ 
summer days would be balanced by short winter days. However, ain^e 
the average zenith angle for the sun is greater at higher latitudes, 
absorption and scattering are more pronounced at such locations. This 
effect alone causes the annual average clear-day values of direct- 
normal irradiance in the northern part of the country to be some 20% 
lower than those in the southern part [A-l]. 
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The geometrical relationship between the sun’s position in the 
sky and the latitude of the observer also leads to greater extremes in 
the number of daily hours of possible sunshine in northern latitudes 
than in southern latitudes. Thus, for example, the extreme lengths of 
a uniat r and winter days at 35°N latitude are approximately 14 and 10 
hours respectively, while at 42°N latitude they are about 15 and 9 
hours? cf, (a) and (b) of Figure A-2. 

D. ATMOSPHERIC TURBIDITY 

Absorption and scattering processes in the earth’s atmosphere may 
be classed conveniently into essentially constant effects and time and 
position dependent effects. The former include absorption and 
scattering due to dry air and ozone and the latter, that due to water 
vapor and dust. 

The term turbidity refers to the amount of absorption and 
scattering by aerosols, and includes effects due to dust particles at 
both high and low altitudes. Both water vapor and low-altitude dust 
concentrations show marked annual as well as geographical variations. 

In each case, the concentration is higher in summer than in winter. 
Water vapor concentration is high for the southeastern states and low 
for the non-coastal northwestern, while low-altitude dust 
concentrations tend to be high in the eastern part of the country and 
low in the western [a-1]. 

Superimposed on this is the effect of high-altitude dust 
particles. The primary source of these is volcanism, and their 
concentration has shown considerable variation over the years. 

Variations in lower atmospheric turbiditv cause direct-normal 
irradiance values in the east to he some 20% lower than those in the 
west, while water vapor concentration variations result in variations 
of only a few percent [A-l]. 

Upper atmospheric turbidity variations caused bv volcanic 
activity result in overall average direct-normal irradiance variations 
of approximately 5-10% with some changes amounting to greater than 30% 
(Krakatoa! [a-1j. 

E. CLOUD COVER 

Cloud formation and distribution show strong temporal and 
geographical variations, and clouds result in a major reduction of the 
direct-normal component. In the absence of clouds, the cumulative 
effect of the phenomena discussed above is a reduction of the average 
clear-day direct-normal component bv about 25-30% in going from the 
southwest to the northeast. The lower percentage of sunshine due to 
clouds results in insolation values in the northeast that are less than 
half those in the southwest and drastically alters the geographical and 
temporal distributions of insolation over the United States. 

Variations in average cloud cover at various locations throughout 
the country show both short-term (monthly! and long-term features. 
Typically, the former amount to some ^ 1 S% while variations over 
periods of a few years are about 10% [a-1]. 
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Figure A-2. Typical Temporal Variation of Solar Insolation on 
Selected Clear Days at Two Different Sites 
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F. 


DISCUSSION 



Owing to the nonstationary nature of the statistics inherent in 
the I DN versus time curve, it is not possible to determine true 
averages for any finite time interval because such average values will 
depend in a random way upon the time selected. Accordingly, caution 
must be exercised in using data based upon averages formed during a 
particular time interval because they may not be representative of 
those formed during a different interval. 

One may nevertheless consider the 1^ curve to have certain 
characteristic time dependences resulting from the phenomena described 
in the previous sections, and an understanding of these is essential 
in considering the influence of solar insolation on solar thermal 
power plant operation. 

Figure A- 3 shows a typical I DN curve for thirty days of one 
month. The non-stationary nature of the statistics is evident in this 
figure. 


IV. INSOLATION DATA 
A. INTRODUCTION 

Although many insolation measurements have been made at a large 
number of sites across the United States and over a period of many 
years, very few data of the kind and accuracy required for the 
engineering analysis of solar thermal power plants are available. The 
main reason for this circumstance is that very few measurements of 
IqN exist; in the past, agricultural and other land use considera- 
tions favored the measurement of 1^ rather than I DN . 

This situation has been further compounded by the fact that large 
errors have been found to exist in much of the total horizontal insola- 
tion data [A-2] and some of the limited direct-normal data available 

0 - 3 ] • 



The following sections summarize the current situation with 
respect to these data. 

B. DIRECT NORMAL INSOLATION DATA 

Since high-efficiency operation of a solar thermal plant dictates 
the achievement of high receiver temperatures with the concomitant need 
for precise focusing of the direct-normal term, there exists a need 
for reliable I ^ data at many locations across the country. 

The measured data which exist may be divided into those which 
have been gathered by the National Weather Service and archived by the 
National Climatic Center (NOAA data*) and those gathered by other 
sources . 


*NOAA - National Oceanic and Atmospheric Administration, Dept, of 
Commerce 
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Figure A-3. Direct Normal Solar Insolations as a Function of 
Time During January 1962 at Albuquerque, NM 
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In addition to the NOAA network stations, there are some 80 
locations where total horizontal measurements are being taken by 
various governmental, state, and private organizations. Some of these 
data are being archived by the NCC, but most are not, and individual 
sources must be contacted in order to obtain information [A-4] . 


V. INSOLATION MODELS 

A* INTRODUCTION 

The paucity of adequate high-quality data has prompted a number 
of workers to develop insolation models based upon different 
assumptions and different kinds of data. In spite of the complexity 
of the problem and the limited measured data available, certain of 
these modeling techniques are quite successful and, in fact, provide 
what is probably the >est means presently available for determining 
accurate insolation v >s on an hourly basis at various locations 
throughout the country. An indication of the success of these modeling 
methods is the fact that one of them is currently being used by the 
NWS in its quality control program for the NOAA network data. Thus, 
if measured data are found to differ from the model predictions by more 
than 5%, the data are re-examined for errors [A-5]. 

The following paragraphs give brief descriptions of the bases 
for several of these models and assess their usefulness with respect 
to solar thermal power plant design and economic studies. 

B. WATT MODEL 

Watt Engineering, Ltd., has constructed an insolation model for 
the direct, total, and diffuse components which is relatively simple 
and quite accurate [A-6] . It is based upon an analysis of the 
physical processes involved in the scattering and absorption of 
radiation by the earth* s atmosphere and takes into account all of the 
phenomena felt to be significant in determining the actual values of 
insolation received at a given location at a given time. 

The model relies upon a minimum of assumptions and, aside from 
basic data on dry-air scattering and absorption, exoatmospher ic 
radiance, and geometrical factors, utilizes as input only data on 
upper and lower atmospheric turbidity, water vapor, and cloud cover. 

One of the strongest features of the model, aside from its 
excellent agreement with measured average insolation values, is that 
it can be used to generate accurate clear-day Ip^ versus time curves 
for the analysis of solar thermal systems. With the addition of a 
more refined handling of cloud statistics, it could also produce such 
curves for cloudy days. Such an approach would enable the construc- 
tion of a realistic data input for any location in the country at any 
given time and, hence, could provide the best means of system 
evaluation for different sites. 



A- 13 



r 




Another strong point is that it does not rely on total insolation 
data of questionable accuracy in order to compute the direct normal 
component. It does, however, rely upon monthly averages for water 
vapor and turbidity data and does not therefore take into account more 
rapid fluctuations in these quantities. 

Figures A-4, A-5, and A-6 present respectively the resu’ts 
obtained by Watt Engineering, Ltd. for direct normal, diffuse 
horizontal, and total horizontal insolation. These figures have been 
reproduced from Ref. A-l. A comparison of the annual average daily 
direct normal insolation characteristics tabulated as Table 1 in the 
body of this report (and reproduced here as Table A-l) with Figure A-4 
shows a fairly close agreement of corresponding values. 

C. AEROSPACE MODEL 

Randall and Whitson [a- 6] have constructed a model for 
determining hourly direct normal insolation values from hourly data 
regarding total horizontal insolation. While the procedure works 
quite well, it suffers from the generally poor quality of the 
available data used as the input. 

The procedure is based upon a statistical analysis of the 
correlation between hourly direct normal measurements and 
simultaneously measured values of total horizontal insolation. The 
data used to construct the model enable the calculation of annual 
average direct normal insolation values from total insolation data 
with an accuracy of 4%, provided the input data are at least this good. 

However, when applied to the rehabilitated NOAA data for 24 
locations in the U.S., the uncertainty in the results is believed by 
the authors to be 10%. Like the Watt model described above, the 
Aerospace model permits the calculation of statistically correct 1^ 
versus time curves for both clear and cloudy days and, hence, may also 
be used for solar thermal system analysis. It is these computed 
hourly values of I DN which appear on SOLMET tapes produced since 
April 1978 for the six NWS stations listed in Table A-l. 

Figure A-7 presents a comparison of annual average daily direct 
insolation contours as determined by the Aerospace Model with those 
determined by the Watt Model. This figure vividly illustrates the 
magnitude of the currently existing uncertainties in an input 
parameter that is critically important to proper site selection and 
subsequent evaluation of the solar thermal power system whose 
operation is estimated at that site. The existence of geographic 
subregions in which annual average daily direct insolation levels are 
especially favorable or unfavorable is suggested by Figure A-7, but 
the figure clearly shows that much better definition of these 
subregions is urgently needed. 
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able A-l. Insolation Characteristics of Chosen Sites 





D. LIU AND JORDAN MODEL 

In the Liu and Jordan approach, observed statistical correlations 
between monthly average values of the daily diffuse horizontal and 
daily total horizontal insolation values, between hourly diffuse 
radiation and solar time, and between hourly total radiation and solar 
time, are used to calculate average hourly values of direct, diffuse, 
and total insolation on surfaces of arbitrary orientation* The input 
data to the model are observed monthly average daily total horizontal 
insolation values for the specific location considered [A-7] . 

While the model yields hourly data, these data are for average 
days and, hence, do not reproduce the temporal effects due to clouds* 
Also, it is subject to the same limitations as the Aerospace model 
through its dependence upon observed 1^ values* Thus, while the 
hourly time dependence is not realistic enough for system dynamic 
analysis, the data may be averaged to yield insolation maps for 
various time periods. 

E. THE HOYT MODEL 

Hoyt [A-8] has recently published model calculations similar to 
but more detailed than those of Watt. These calculations use measure- 
ments of turbidity, precipitable water vapor, and surface albedo of 
the earth as input data, and result in clear day solar noon total 
hemispheric insolation values that agree with calibrated pyranometer 
readings to within 2. 72. The agreement between the model and the 
measurements is well within the instrumental error of ±5%. 

These calculations have been used by NOAA to rehabilitate the 
total insolation data from 26 NWS stations in the network and, as 
mentioned in the introduction to this section, are also used as a 
means of quality control of the data. 

The model has also been used to calculate daily average 
clear-day values of as well as monthly mean values of I^h 
( including clouds). In all cases the agreement with measured values 
is within the pyranometer error. 

While the model has good accuracy, it does not permit the calcu- 
lation of hourly values of 1^, either with or without clouds, nor 
does it permit the calculation of Iq^, although it could be modified 
in such a way as to include these cases. 


F. METHODS FOR GENERATING CONTOUR MAPS 

The question of the effect which different interpolation schemes 
have on contours produced from a given data set does no" appear to 
have been considered by those presenting such maps because little or 
no mention is generally made as to which of the many possible methods 
for generating contours has actually been used. The discussion which 
follows presents some of the basic problems that should be considered 
and provides a number of specific examples which demonstrate the 


A- 20 


r 


rather large variations that can be obtained by using different 
methods for the generation of contour maps. At present, there does 
not appear to be any generally accepted criterion which would clearly 
indicate which method should be chosen to the exclusion of all others. 

Two basic approaches may be used to generate contours from a 
given set of randomly distributed points at which z-values are 
specified. In the first, the points are interconnected to form a net 
of triangles, and points on the contour lines are determined by linear 
interpolation of the z-values along each of the sides of each 
triangle. Common z-value points are then joined by a smooth curve. 

The second approach begins by fitting a smooth surface to the z-values 
given, and then drawing the contours. 

In the first approach, there is ambiguity in the choice of 
interconnections in that there is no unique way to join a set of 
randomly distributed points to form triangles. Thus, different sets 
of triangles will result in different contours. Also, any of several 
different methods can be used to draw smooth curves through the common 
z-value points found by interpolation and, generally speaking, each 
different method will result in a slightly different set of contours. 

In the surface fitting approach, all of the ambiguity is 
concentrated in the surface fitting algorithm since, once a smooth 
surface is obtained, the contours are unique.* 

The number of wavs a smooth surface can be fitted to a randomly 
distributed set of z-values is virtually unlimited. Some common 
methods employ splined cubic functions while others adapt an iterative 
approach in which the surface is constrained by some simple partial 
differential equation (Poisson's equation, the biharmonic equation^. 
Even with a given method, the resulting surface and, hence, the 
contours depend upon the exact manner in which the algorithm is 
applied. For example, one method of fitting a splined cubic function 
involves the initial formation of a triangular net, as in the linear 
interpolation scheme described earlier. Since the polygon formed in 
the (x,y) plane must be convex in this method, one must usually supply 
additional data points to accomplish this, and since these points must 
themselves be determined by some method of interpolation or 
extrapolation, an additional uncertainty is introduced into the final 
contour map. Actual computer runs show that such schemes tend to 
produce questionable edge effects which depend upon the exact choice 
made for assigning z-values to the extra points. 

Another important consideration that enters into the surface- 
fitting problem is the extent to which the final surface actually 
coincides with the original data points. Many techniques of smoothing 
introduce departures from a perfect fit unless the algorithm 


♦Actually, one may view the triangle method as a crude surface fit 
which would produce a unique set of linearly segmented contours. 
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explicitly prevents this by constraining the surface to pass through 
the original points as the iteration proceeds. Such departures nay or 
may not be desirable depending on the data and the purpose for which 
the contours are to be used. For example, if, as is often the case, 
the existing data constitute a severely under-sampled set, a smoothed 
surface which does not necessarily provide a perfect fit might be more 
representative of average trends than one which does. For example, a 
random sampling that happens to select adjacent high and low z-values 
would lead to unrepresentative contours if the surface were 
constrained to pass through these points. 


Tht! above discussion suggests an approach to surface fitting 
based upon Fourier theory analogous to that used when the data points 
form a periodic lattice. Such an approach is briefly described in the 
following paragraphs. 


A given set of randomly distributed z-values may be represented 
by the following sum, 


where: 


N 

s(x,y) ■ £ *£ Mx “ O My - yM, 

>1 

z; * value of Z(x,y) at the point (x{, yj), 

6(x) ■ Dirac delta function, 

N ■ number of data points. 


The Fourier spectrum of this set of weighted 6-functions is: 
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where: 

K x “ spatial frequency in x-direction, 

► y - spatial frequency in y-direc t ion. 

If the average separation between data points in the x-direction 
is Ax, and that in the y-direction is Ay, and we choose to construct a 
surface z'(x,y) such that the given points constitute an average 
Nyquist sampling of this surace, then the maximum spatial frequencies 
that the Fourier spectrum of the surface can have in the x and y 
directions are: 


k 

xm 


Ax 
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Thus, if we apply a low pats filter, F (► x , ► , ) , having cutoff 
frequencies of and ► , n to the apectrua of our original aet of 
weighted 6-function* , we obtain: 

r'(^ , ► ) - ?(< , ► ) ?(< , O , 
x’ y x’ y x' y ’ 

and taking Fourier transform* we find: 


r 

l 


where: 


N 

r'(x, y) - ZjF(x - x £ , y - y £ ) , 

i"l 


/V 

F(x, y) * Fourier transform of F(* x , ► y ). 

We are thua left with a surface formed by replacing each original 
x-value by the spatial .1 filter F(x,y) multiplied by that 

z-value. For example, if F (x x , k y ) is a Gaussian function, then 
F(x,y) is also a Gt,uS»ian, and the surface z'(x,y) is a weighted sum 
of Gaussians. Cleuriy, our filter should be chosen so that these 
overlap as little as possible if the surface is to pass through the 
original z-values, and we must also choose our filter such that 
F(0,0) * 1, that is, we must demand that 


and 



► ) d> d 

x y 


1 


F(x n - x{, y n - y{) *> 1 $ for all i, n. 

The above procedure can be expected to provide a representative 
surface fit provided there is r.o clustering of points in the (x,y) - 
plane. A clustering will cause strong overlap of the functions 
F(x - x^, y - y{), resulting in abnormally high values in the 
region where the clustering exists. Also, the choice of a filter 
function must be made so as to avoid ringing on the one hand 
(F(* x , -y) t0 ° steep near cutoff) and a non-smooth surface on the 
other (Ftx,y) a box filter, for example). A Gaussian function appears 
to be a natural choice. The clustering problem can be handled by 
departing from linear filter theory and directly applying a spatially 
dependent filter in the spatial domain by means of convolution. Then, 
each point will have its own filter chosen so as not to overlap those 
of neighboring points, thus giving explicit recognition to the fact 
that the sampling represented by the original data set is non-uniform, 
i.e., some regions contain more information than others. 
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In order to assess the confidence that one should have in an 
insolation contour map produced from a given set of data by a given 
surface fitting/contouring technique, a series of maps has been 
prepared from two different data sets by a number of different 
techniques so that the reader may compare these and judge for himself 
what limits to place on conclusions drawn from one map or another. 

The two data sets used in the study are the Watt Engineering 
Ltd. estimates of direct normal insolation, and the Aerospace Corp. 
estimates of direct normal insolation, cf., Sections B and C above. 
Annual averages based on approximately 25 years of recent data have 
been used in producing both the Watt and Aerospace contour maps. 
Figures A-8 and A-9 show a comparison of contours produced from Watt 
and Aerospace data by means of simple linear interpolation. Although 
the general trends are similar, the two maps differ considerably in 
detail. Figures A- 10 and A- 11 show contours produced from the same 
two data sets by fitting a splined cubic function to the given data 
points. Again, there is considerable difference in the contours. A 
comparison of Figure A-8 with A-10, and Figure A-9 with A-ll shows the 
different effects which these two contouring methods have on the 
results, and, while the differences are smaller than those due to the 
different data sets, they are still quite apparent. 

Figures A-12, A-13, and A-14 present the results of applying 
various low pass filters to surfaces produced from the Watt data by a 
growth procedure wherein a discontinuous surface is first produced by 
assigning to each (x, y) point the z value corresponding to the 
nearest data point. It is seen here that one may obtain significantly 
different contours depending upon the degree of smoothing provided by 
the filter. The intermediate filter (Fig. A-13) corresponds 
approximately to an average Nyquist sampling, and the problem 
associated with clustering has been handled by applying low pass 
filters to the discontinuous surface produced by the growth procedure 
rather than to the set of weighted 8 -functions represented by the 
original data. This prevents excessively high values from appearing 
in the final surface and permits the use of a simple box filter. 

Figures A-15, A-16, and A-17 show the results of applying the 
same technique to the Aerospace data, where Figure A-16 corresponds to 
an average Nyquist sampling. The effect of the higher sampling rate 
of the Watt data is clear from a comparison of Figures A-12, A-13, and 
A-14 with Figures A-15, A-16, and A-17. 

The maps shown in Figure A-12 to A-17 were generated at the 
Image Processing Lab (IPL) at JPL, and illustrate the power of thi9 
facility in producing such contoured surfaces. An additional benefit 
provided by IPL techniques can be seen in the continuous gray scale 
representation shown in the figures; this gray scale provides a 
graphic display of the surface, and, hence, of the actual insolation 
values represented by the filtered data. Such contour maps can be 
generated quickly and inexpensively once the original U.S* map has 
been produced and stored; thus, new data may be incorporated readily 
into new maps. 
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Figure A-9. Annual Mean Daily Direct Normal Insolation 
(Aerospace data, linearly interpolated) 
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i Daily Direct Normal Insolation 
data, cubic spline fitted) 
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Figure A-i A . Data Based on Watt Engineering, Ltd. Estimate of Direct Normal Insolation. The Degree of 
Smoothing of the Surface Increases in Going from Figure A-12 to Figure A-14, Figure A-13 
Corresponding Approximately to an Average Nyquist Sampling (See Text). 
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of the Surface Increases in Going from Figure A-15 to Figure A-17. Figure A-16 Corresponds 
Approximately to an Average Nyquist Sampling (See Text). 





Figure A-17. Data Based on Aerospace Corp. Estimates of Direct Normal Insolation. The Degree of Smoothing 
of the Surface Increases in Going from Figure A-15 to Figure A-17. Figure A-16 Corresponds 
Approximately to an Average Nyquist Sampling (See Text). 
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DISCUSSION 


The above is not an exhaustive discussion of existing insolation 
models; nevertheless, it provides an indication of the range and 
accuracy of the calculations that have been made and indicates which 
approaches are most useful for solar thermal power systems analysis. 

In this connection, one must consider two distinct needs 
requiring two essentially different kinds of data. On the one hand, 
systems analysis, either engineering or economic, requires at least 
hourly data which simulate the actual insolation values (Idn or 
I?h) with correct statistics for given locations. On the other 
hand, site selection considerations require the construction of 
insolation maps based upon longer term averages (monthly, annual), 
with provision for estimating long-term variations. Thus, somewhat 
different approaches must be used in each rase. 

Figure A-18 is included to show the very different distribution 
in monthly average daily direct insolation that may exist at one site 
as compared with another, even though the annual average daily direct 
insolations for the two sites may be nearly the same. At each site, a 
particular recent year was selected on the bases that (a) the data 
were recent enough to insure reasonably accurate measurement and (b) 
the annual and monthly average daily direct insolation values were 
close to the corresponding long-term averages. Referring to Figure 
A-18, the distributions of monthly average daily direct insolation at 
Medford, OR, and Miami, FL are drastically different. During June, 
July, and August, average daily direct insolation at Medford, OR, is 
virtually equal to that at Albuquerque, NM; during November, December, 
and January, average daily direct insolation at Medford, OR, (Ij)N « 

.9 kWh/day) is several times lower than that at any of the other sites 
for which data are presented here. On the other hand the monthly 
average daily direct insolation is virtually constant throughout the 
year at Miami, FL. 

Figure A-19 presents, for each of the sites shown in Figure 
A-18, the number of hours during the selected year that hourly average 
direct insolation was between specified intensity limits defined by 
increments of .05 kW/m^. These histograms complement the results 
presented in Figure A-18. The relatively even distributions of 
intensity at Miami and Maynard reflect the climatological character of 
the Eastern coastal climate. Madison, as well as Miami and Maynard, 
shows a sharp decrease in incidence of the higher intensity levels 
(0.85 - Idn^J this attributed to the high relative humidity 
established at these sites by surrounding bodies of water and frequent 
rainfall whenever the local weather is warm. On the other hand, the 
histograms for Barstow and Albuquerque reflect the large proportion of 
days that are clear and during which, because the surrounding regions 
are relatively arid, the relative humidity is low. These 
characteristics, coupled with low latitudes, appear to be the chief 
factors contributing to the relatively much higher performance of 
solar thermal power systems located in the Southwest as compared with 
these same systems located in other regions of the contiguous United 
States . 
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Figure A-18. Monthly Average Daily Direct Insolation lor a 
Typical Year at Each of Several Sites Within 
the United States 
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r igure A- 19. Histograms of the Annual Hourly Average hi root 
insolation lor the Same Years and tor the Sana 
Sites Presented in Figure A- 1 8 . 
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Figure A-19 (cont'd) Histograms of the Annual Hourly Average Direct 

Insolation for the Same Years and for the Same 
Sites Presented in Figure A-I8. 
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SUMMARY 


The basic physical considerations relevant to insolation have 
been presented, including a discussion of the quantities measured and 
the basic instrumentation used for such measurements* 

Various factors influencing the temporal and geographical 
variation of insolation values have been analyzed, and estimates of 
the variations caused by these factors have been given* Across the 
contiguous United States, these include a north-south variation of 
20% due to zenith angle effects, an east-west variation of 20% due to 
atmospheric turbidity variations, and a 2-to-l southwest-northwest 
variation in average cloud cover. In addition, the fact that the 
dominant long-term effect of upper atmospheric turbidity caused by 
volcanic activity can result in changes in northern hemispheric 
insolation values of more than 30% within a year or so, has been 
emphasized. 

According to recent studies (Ref. A-9), insolation measurements 
have nominal error bands of about ±3%. For coastal regions at 
extreme northern latitudes (e.g., Seattle, Washington, and Caribou, 
Maine), it appears that actual direct insolation could be of the order 
of 20% higher than is indicated by SOLMET insolation tapes. This 
would tend to improve performance of solar plants in these extreme 
northern coastal regions but would not materially affect the 
comparisons among plant types. Moreover, performance in these 
northern, low-insolation regions would still be lower than that in the 
Southwest and other lower-latitude regions. 

The existing insolation data base in the United States has been 
analyzed with respect to accuracy, completeness, and type of data, and 
it has been pointed out that this data base is not completely adequate 
in each of these respects for every solar thermal study. 

The current status of model calculations has been investigated 
and the conclusion reached that certain of these models, individually 
or in combination, are adequate for calculating any type of insolation 
at any location and for any given time frame. Thus, they are suitable 
for both quasi-steady-state systems analysis and site selection 
studies . 


VII. inclusions 

From the above it is clear that, in the future, the evaluation 
of solar thermal systems must depend upon accurate model 
calculations. None of the existing models are entirely satisfactory 
for this task at the present time, but, taken as a whole, they provide 
a basis upon which a suitable model can be constructed. A more nearly 
ideal model would 

1. permit the calculation of I DN , I TH , and values; 

2. allow such calculations to be made at any specific site 
within the continental United States; 
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3. allow such calculations to be made for any given time 

period and to any practical time resolution (time interval 
for integration) ; 


4. do all of the above with correct statistics. 

In order to construct such a model, the following steps need to 
be taken: 

1. Improved and more extensive measured data must be acquired. 

2. The existing models must be analyzed and compared for 
overall suitability in terms of accuracy and flexibility. 

3. One or more of these must be selected as a basis for 
further modification to allow the proper inclusion of the 
dynamic effects of clouds. 

4. The resultant model must be checked for accuracy and 
evaluated with respect to sensitivity to input data. This 
would involve surveying the existing sources of such data 
as turbidity, water vapor, and albedo and determining the 
errors in insolation values introduced by errors in these 
parameters. 


A-40 


REFERENCES 


A-l. "On the Nature and Distribution of Solar Radiation," DOE Report 
HCP/T2552-01. Watt Engineering Ltd-, (Mar. 1978). 

A-2. "Solar Radiation Data Sources, Applications and Network Design," 
DOE Report HCP/T5362-01. University of Alabama, pp. 1-18, (Apr. 
1978). 

A-3. Robert W. Durrenberger and Anthony J. Brazel, "Need for a Better 
Solar Radiation Data Base,” Science lj), 1154-5, (Sept. 1976). 

A-4. Carter, E. A., Wells, R. E., and Williams, B. B., "Solar 

Radiation Observation Stations with Complete Listing of Data 
Archived by the National Climatic Center, Asheville, N.C., and 
Initial Listing of Data not Currently Archived," Center for 
Environmental and Energy Studies, the University of Alabama, 
Huntsville, Nov. 1976 - prepared for ERDA, Div. of Solar Energy, 
under contract NA58-31293. 

A-5. Ward Seguin, National Weather Service, (private communication). 

A-6. "A New Procedure for Estimating Hourly Direct-Normal Insolation 
from Hourly Total-Hemispheric Insolation," C. M. Randall and 
M. E. Whitson, Aerospace Report ATR-78(7592)-l, Dec. 1, 1977. 

A-7. "Hourly Solar Radiation Dat.. for Vertical and Horizontal 

Surfaces on Average Days in the United States," T. Kusuda and 
K. Ishii, NBS Building Science Series 96, U.S. Dept, of 
Commerce, (April, 1977). 

A-8. Douglas V. Hoyt, "A Model for the Calculation of Solar Global 
Insolation," Solar Energy 21, 27-35, (1978). 

A-9. Rapp, D., "A Critique on Rehabilitation Procedures Used by NOAA 
for Solar Data," to be published in Energy Conversion Journal 
(1979). 


A-41 


APPENDIX B 

PERFORMANCE AND COST CHARACTERISTIC! 
OF 

ADVANCED SOLAR THERMAL POWER PLANTS 


B-l 


CONTENTS 


I. INTRODUCTION B-4 

II. SYSTEM CHARACTERIZATION B-4 

A. PERFORMANCE DATA B-A 

B. COST DATA B-8 

C. FINANCIAL FACTORS B-12 

III. OPERATION IN REGIONS WITH DIFFERING INSOLATION B-13 

A. COMPUTER SIMULATION METHODOLOGY B-13 

B. MINIMUM ENERGY COST ENVELOPES B-1A 

C. COLLECTOR FIELD AREAS AND STORAGE TIMES B-19 

D. CAPITAL COST BREAKDOWN FOR SELECTED SYSTEMS B-26 

IV. REFERENCES B-29 


IV. REFERENCES B-29 



FIGURES 


B-l . System Design Efficiency Train B-5 

B-2. Part Load Engine Efficiencies — - — - — B-9 

B-3. Structure of Power Plant Simulation Computer 

Program — B-l 3 

B-4. Minimum Energy Cost Envelopes for the Paraboloidal 

Dish Stirling System B-15 

B-5. Minimum Energy Cost Envelopes for the Central 

Receiver System B-16 

B-6. Minimum Energy Cost Envelopes for the Parabolic 

Trough System B-l 7 

B-7. Effect of Regional Insolation Differences on 

Collector Field Area for Capacity Factor • 0.55 B-20 

B-8. Effect of Regional Insolation Differences on 

Collector Field Area for Capacity Factor * 0.40 B-21 

B-9. Effect of Regional Insolation Differences on 

Collector Field Area for No Storage B-22 

B-10. Effect of Regional Insolation Differences on 

Storage Time for Capacity Factor ■ 0.55 B-23 

B-ll. Effect of Regional Insolation Differences on 

Storage Time for Capacity Factor ■ 0.4 B-24 

B-12. Effect of Regional Insolation Differences on 

Capacity Factor for No Storage B-25 

B-13. Capital Costs for Selected Power Plants at 

Capacity Factor “ 0.55 B-27 


B-3 


APPENDIX B 


I . INTRODUCTION 


This appendix first describes the rationale and procedures used 
to determine the performance and cost characteristics of the four 
selected advanced solar thermal power plants. 

o Paraboloidal Dish (PD) 

o Central Receiver (CR) 

o Parabolic Trough (PT) 

o Compound Parabolic Concentrator (CPC) 

Then, selected results based on computer simulations of these systems 
operating in regions with differing insolation patterns (Appendix A) 
are presented. These results provide details which supplement and 
amplify the primary findings summarized in the body of this report. 


II. SYSTEM CHARACTERIZATION 


The four selected systems were characterized in terms of 
performance and cost data pertaining to the subsystems/components 
comprising each system. The ground rules used in determining 
performance and cost characteristics were: 

o Subsystem/component advanced technology development was 
projected to the 1990-2000 timeframe. 

o High-volume mass production ( 10^ units/year) was 

assumed for modular components (e.g., concentrators). 

These ground rules are closely related to those used in 
determining solar thermal program targets. Thus, a basic source of 
information regarding the desired projected data is program planning 
documents (e.g., Refs. B-l and B-2). These planning documents were 
used as the starting point. In general, the information available was 
not adequate for the detailed system characterization required for 
computer simulation methods employed in the study. Therefore, it was 
necessary to generate much of the detailed subsystem/component 
performance and cost data base. The basic approach was to work with 
specialists from Government Laboratories and industry who are actively 
involved in developing the four selected systems and their components. 


A. PERFORMANCE DATA 

The performance of the selected power systems is presented in 
Figure B-l. It shows the design-point efficiencies of components in 
the primary power train that captures solar energy and converts it to 
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Figure B-l. System Design Efficiency Train 


electricity. This electricity is then delivered to the grid or 
application site. It is shown that the two-axis tracking systems 
represented by the parabolodial dish and central receiver achieve 
higher design point system efficiencies than the one-axis parabolic 
trough. The paraboloidal trough has a higher efficiency than fixed or 
periodically adjusted orientation systems such as the Compound 
Parabolic Concentrator. 

These basic trends explain the motivation behind efforts to 
develop tracking concentrator systems. Tracking permits a greater 
fraction of the available insolation (solar flux) to be captured. 

This solar flux can be highly concentrated at the receiver, 
particularly in the case of two-axis tracking systems. The result is 
higher temperatures and correspondingly higher conversion efficiencies 
as seen in Figure B-l. 

For all systems, the reflective surface was taken to be a 
back-silvered, thin, low-iron content glass produced, e.g., by the 
fused glass process. Based on laboratory tests at Sandia 
Laboratories, Albuquerque, (Refs. B-3 through B-6) reflectivities in 
the range of 0.90 to 0.95 are projected. 

In addition to reflectivities, each of the systems requires the 
specification of key design parameters. Slope errors for the 
paraboloidal dish of ±0.1° are assumed for the mirrored reflective 
surface* Low slope errors yield high concentration ratios with small 
spillage of flux outside the receiver aperture. The result is 
efficient capture of solar flux for operation at high temperatures in 
the range of 1500°F to 2000°F. These temperatures are deemed 
appropriate for power generation from paraboloidal dish systems (Ref, 
B-7 ) . 


A major consideration for the central receiver is the layout of 
the heliostat field with respect to the tower. The highest optical 
efficiencies for the latitudes encompassing the contiguous United 
States (Refs. B-8 and B-9) are provided, by a north-field arrangement, 
wherein the heliostat field is located north of the tower. However, 
the north-field requires a higher tower for the receiver than do field 
arrangements where the tower is more centrally located within the 
field. As the power rating and size of the field are reduced, the 
tower height decreases. Tower design complexity and costs then 
diminish and north-field designs appear to be preferable for small 
central receiver systems. In the present study, an efficient 
north-field arrangement was assumed even for a 10 MW plant. The 
higher cost of the north-field tower was offset by projecting 
development of advanced high-temperature liquid metal or molten-salt 
heat transport systems to carry energy from the tower-mounted receiver 
to the energy conversion system. The heat transfer characteristics of 
liquid metals like sodium allow construction of lighter weight 
receivers and the use of smaller pipe diameters as compared to steam 
or gas systems. This results in the projections of lower tower costs 
according to estimates by Sandia Laboratories, Livermore (Ref. B-10). 

Minimization of receiver heat losses is required for high 
efficiency operation of the parabolic trough. For a reflective 
surface with a small slope error (10.1°), the reflected solar flux 
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will be concentrated in a narrow strip impinging on the linear 
receiver. To minimize heat losses, it is necessary to insulate the 
receiver; on the other hand, a thick insulation will block a portion 
of the insolation from the concentrating surface. Thus, sophisticated 
insulation techniques as, for example, evacuated jackets in 
conjunction with advanced insulation materials will be required. 
Performance projections for the trough are predicated on evolutionary 
development of this technology as supported by activities at Sandia 
Laboratories, Albuquerque. These projections indicate that the next 
generation of parabolic trough systems will perform much more 
efficiently than do present systems (Ref. B-ll). 

The CPC system (Refs. B-12 and B-13) is less critically 
dependent on reflective surface accuracies and can capture a portion 
of the diffuse radiation. A key parameter for this system is the 
concentration ratio. As the concentration ratio is increased by 
enlarging the reflecting surfaces, higher operating temperatures are 
practicable. However, as the sun passes over the fixed reflecting 
surfaces of the CPC, only the direct insolation falling within the 
acceptance angle or field-of-view is concentrated on the receiver. 

The acceptance angle becomes smaller as the concentration ratio is 
increased. Based on guidance from Argonne National Laboratory, a 
concentration ratio of 3 to 5 was employed in this study. Sensitivity 
to concentration ratio, particularly for CPC systems operating in 
regions with high percentages of diffuse radiation, is discussed 
further in Appendix C. 

The efficiency of transporting energy is high (>90%) for all 
systems, even though widely different approaches are represented. The 
modular paraboloidal dish relies primarily on electric transport, 
whereas the central receiver employs transfer of energy from 

the field. The parabolic trough and CPC use thermal transport of 
steam and organic fluids, respectively, via pipelines. The central 
receiver, with an optical atmospheric attenuation loss of only ~2% for 
a 10 MWe plant, has the highest projected transport efficiency. The 
assumed values for electrical and thermal transport efficiencies could 
theoretically be improved, e.g., by using larger diameter wires to 
reduce electrical resistance losses and larger diameter pipes with 
thicker insulation to reduce frictional and thermal losses. However, 
these measures would increase costs. Efficiency values are therefore 
a function of engineering trade-offs. For this study, values 
consistent with current electrijal and pipeline networks were u^ed 
since they reflect engineering trade-off considerations. 

A spectrum of energy conversion systems is employed in the 
selected advanced systems. For the paraboloidal dish, a Stirling 
engine operating at 1650°F is employed. For the central receiver, 
an 1800°F combined-cycle Bray ton/Rankine system is used. The 
parabolic trough employs a Rankine cycle steam power system operating 
with steam at 800°F and the CPC employs an organic Rankine system. 
Current engines generally achieve about 60 percent of the 
corresponding Carnot efficiency. Carnot efficiency is the maximum 
achievable for specified conditions. With the development of advanced 
designs and the incorporation of features such as reheat and 
intercooling, it is possible to improve the efficiency of all the 
engine types above. 


Trade-offs between improved efficiency, added cost, and 
complexity will dictate the extent to which efficiency is improved. 

For solar thermal plants, the achievement of high efficiency is 
extremely attractive since the collector field area can be reduced as 
the efficiency increases. System trade-offs favor development of 
efficient engines and, in this context, it is projected that values of 
about 70 percent of Carnot cycle efficiency will be attained in the 
1990-2000 timeframe. Since Carnot efficiency increases with operating 
source temperature, the higher source temoerature systems have higher 
conversion efficiencies, as seen in Figure B-l. 

A value of about 70 percent Carnot cycle efficiency is the 
practical maximum design efficiency for the energy conversion systems. 
The basic types of engines have differing part-load efficiencies as 
shown in Figure B-2. All the systems have excellent part-load 
performance in that they maintain high efficiencies over a wide load- 
range (Ref. B-7). The Bravton/Rankine achieves its peak efficiency at 
a power less than its design power rating. It is important to 
understand that, only at its maximum normalized efficiency of about 
1.02 (occurring at approximately 0.7 of design power output), does the 
Brayton/Rankine power conversion system associated with the CR achieve 
its practical design maximum corresponding to 70Z of Carnot cycle 
efficiencv. This curve is characteristic of closed-cycle, 
recuperated, Bravton engines (Ref. B-14). 


B. COST DATA 

Each of the four selected advanced power plants was separated 
into the following cost elements: 

• Land 

• Concentrators 

• Receivers 

• Energy Transport 

• Storage Structure/Med ium 

• Heat Engine Equipment 

• Balance of Plant 

• Operations and Maintenance 

Unit cost data for each of these elements was determined bv (1) using 
cost target data from planning documents, (2) synthesizing results 
from mass-production studies, and (3) conferring with personnel from 
other Government Laboratories who are responsible for developing the 
selected systems. Data for the paraboloidal dish Stirling (PDS) 
system were determined within the Jet Propulsion Laboratory ( JPL) . 

For the central receiver (CR) data inputs, Sandia Laboratories, 
Livermore, the Electric Power Research Institute (EPRI), McDonnell 
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NORMAIIZED EFFICIENCY 



Figure B-2. Part Load Kngine FHicieneios 






Douglas Astronautics, Inc* (MDAC), the University of Houston, and the 
NASA Levis Research Center (LeRC) provided valuable assistance to 
JPL« Similarly, Sandia Laboratories, Albuquerque (SLA) and the 
Argonne National Laboratory (ANL) gave helpful assistance to JPL in 
establishing the data for the parabolic trough and the compound 
parabolic collector (CPC) respectively. 

With regard to land costs, it is tacitly assumed that solar 
plants will not be located on premium real estate since large areas 
are required. Instead, it is considered likely that solar plants will 
be located on low cost land ranging from $1000/acre to $5000/acre. 

The lover coat of $1000/acre is assumed for the arid sunbelt which 
contains large areas of low-cost desert land. 

High-volume mass production studies have only been conducted for 
the central receiver system (Refs. B-15 and B-16). Comparable studies 
at the same production rates are not available for the other systems. 
Therefore, the central receiver studies were used as a baseline. In 
general, existing cost estimates for the collectors of both the PT and 
CPC under limited mass production were higher than those for the CR. 
However, both of these systems were judged to be simpler to 
manufacture and assemble than is the two-axis tracking central 
receiver. Both SLA and ANL concurred with this judgment. Thus, lower 
collector unit costs were used for the trough and CPC, with the CPC 
having the lowest unit cost since it does not require tracking 
mechanisms. 

The paraboloidal dish system, on the other hand, was judged to 
be more complex than the heliostat (concentrator) of the central 
receiver. The dish requires doubly curved reflective surfaces and 
additional structure to support the receiver/engine-generator assembly 
at its focal point. 

The following nominal collector costs have been used to define 
collector costs as the sum of the concentrator and receiver costs. 


System 

Paraboloidal Dish Stirling 
Central Receiver 
Parabolic Trough 
Compound Parabolic Cone. 


Col lector Cost , $/m 2 


95 

75 

70 


These nominal values are considered to reflect relative cost trends. 
The effect of uncertainties in projecting these nominal values is 
evaluated via sensitivity analyses in Appendix r . 

Tor energy transport, the central receiver employs a relatively 
short length of pipe for transporting thermal energy in the tower, 
whereas the parabolic trough and CPC use extensive piping networks to 
bring thermal energy from the collector field to a centrally located 
energy conversion unit. The paraboloidal dish system uses an 
electrical network to collect the energy generated by the small power 
conversion units that are located at the focal point of each dish 
concentrator. 
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Transport costs for the central receiver are included in the 
receiver/ tower subsystem. Costs for the tnermal piping networks of the 
parabolic trough and CPC were estimated on the basis of piping layouts 
for designs being evaluated and studied by SLA and ANL, respectively. A 
pipe network optimization code developed at JPL was used as an aid in 
selecting nominal costs (Ref. B-17). Activities are underway at JPL to 
update electrical network transport costs for the paraboloidal dish 
system. However, this effort is not yet complete; therefore, values 
from earlier studies (Ref. B-7) have been used wherever these values are 
considered to represent nominal \alues. 

Storage costs were inferred from studies (Refs. B-18 through B-22) 
which investigated storage systems in the context of advanced technology 
developments. The advanced battery costs of $30/kWeh for the PDS 
system is consistent with advancid technology projections in Ref r B-19. 
Thermal storage costs are based on projections for advanced latent heat 
systems. Since the trough and CPC employ lower temperature storage than 
the central receiver, lower storage costs are assumed. Containment 
material problems are generally less at reduced temperature levels. 
Recent thermal storage cost goals based on NASA/LeRC activities yield 
lower cost than those goals used in this study, and, therefore, 

Appendix C treats systems sensitivity to storage cost uncertainties. 

Based on high-volume mass production practices as employed in the 
automotive industry, a nominal target cost for the power conversion 
system has been set at 60 $/kWe (Ref. B-7). Automotive industry 
practice pertains to small engines as used in the paraboloidal dish 
system. For the larger 10 MWe plants used in the other systems, the 
same level of mass production will not be achieved if the same overall 
level of market penetration is assumed for each system. Further, 
procedures for mass-producing large power systems are not well 
established. 

However, in the context of a significant market penetration in the 
1990-2000 timeframe, production of large power systems could benefit to 
some extent from increased production. In view of these possibilities, 
it is believed that, as an ultimate goal, large engine system costs 
could approach those of smaller engines. Therefore, the same cost of 
$60/kWe was used for these engines. The increase in energy cost that 
would result if these ultimate values for large engines are not met is 
presented in Appendix C as part of the sensitivity analysis. 

Balance-of-plant includes items such as buildings, electrical 
plant equipment, and master controls. The cost of all these items will 
depend on detailed plant design considerations. For nominal projection 
purposes, the same value of $85/kWe determined by recent central 
receiver studies has been assumed to hold approximately for all the 
svstems. For svstems with large central engines, the plant and engine 
rating are the same but, for the PDS systems using battery storage, 
tote" system engine rating exceeds plant rating. For the present study, 
the balance-of-plant costs were based on total system engine rating for 
the PDS system. For dish systems having high capacity factors and 
correspondingly large field areas, total system engine rating is 
significantly larger than plant rating. At large capacity factors, this 
results in much larger balance-of-plant costs for PDS systems. This 
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treatment of balance-of-plant costs is considered to be conservative 
for the dish system since some balance-of-plant items, such as 
buildings, probably will not increase in proportion to total system 
engine rating. 

Moreover, one basis of analysis used for the PDS system in this 
report is to assume that any power produced in excess of plant rating 
must be stored or wasted. Thus when storage is filled, this basis 
results in the wastage of electrical energy. On an alternate basis, 
the wastage would be fed to the associated electrical distribution 
network and could therefore be credited to the PDS system. 

The estimation of Operation and Maintenance (O&M) costs is 
presently subject to considerable uncertainty. Based on EPRI studies 
(Ref. B-23), a value of about 1% of capital costs is assumed for plant 
operation and maintenance during the first year of operation. The 
effect of variations in O&M costs is treated in the sensitivity 
analysis of Appendix C. 


C. FINANCIAL FACTORS 

The energy cost for solar power plants is determined by using a 
levelized cost methodology developed for utility systems (Ref. B-24). 
This methodology requires financial inputs. Recently, a set of 
financial factors has been adopted for solar thermal comparison 
studies being performed by the Solar Energy Research Institute (SERI), 
Battelle Pacific Northwest Laboratories (BPNL), and JPL. These 
factors listed have been used in this study. 

Utility Description Data 


System Operating Lifetime 30 years 

Annual "Other Taxes" as a fraction of capital 0.02 

Annual Insurance Premiums as a fraction of capital 0.0025 

Effective Income Tax Rate 0.40 

Ratio of Debt to Total Capitalization 0.50 

Ratio of Common Stock to Total Capitalization 0.40 

Ratio of Preferred Stock to Total Capitalization 0.10 

Annual Rate of Return on Debt 0.08 

Annual Rate of Return on Common Stock 0.12 

Annual Rate of Return on Preferred Stock 0.08 

General Economic Conditions 

Rate of General Inflation 0.06 

Escalation Rate for Capital Costs 0 06 

Escalation Rate for Operating Costs 0.07 

Escalation Rate for Maintenance Costs C 

Base Year for Constant Dollars 1. 

Nominal Intermediate Outputs 

Cost of Capital to (and internal rate of return in) 0.086 

a Utility 

Capital Recovery Factor 0.0°39 

Annualized Fixed Charge Rate 0.1568 
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III. OPERATION IN REGIONS WITH DIFFERING INSOLATION 


Computer simulation results for the selected four advanced so'ar 
thermal power systems operating in regions with differing insolation 
characteristics are determined by using the performance and cost 
values of the previous section as inputs. First, the computer 
methodology of simulating solar thermal power plants is described. 
Then, selected results from the computer simulation are presented. 


A. COMPUTER SIMULATION METHODOLOGY 

The structure of the power plant computer simulation program is 
illustrated in Figure B-3. Detailed documentation for this program is 
given (Ref. B-25). The basic input to the program is hourly 
insolation data which are provided in the form of SOLMET tapes for a 
number of sites within the Continental United States (see Appendix A). 
Plant power rating and performance characteristics of 
subsystems/components are the basic inputs to the plant performance 
simulation subprogram. Performance is specified in terms of the 
design efficiency train of Figure B-l and off-design characteristics, 
e.g., engine part-load efficiency as given in Figure B-2. 

Component cost characteristics comprise the basic input to the 
plant capital costs subprogram. These cost characteristics are 
generally unit costs normalized to concentrator area or power rating 
as presented in the previous section. To determine plant capital 
costs, these unit costs are multiplied by collector area, storage 
capacity, and relevant power levels (thermal or electrical) within the 
plant as determined by the performance subprogram. 


PLANT POWER RATING 

PERFORMANCE 

CHARACTERISTICS 


COMPONENT COST 
CHARACTERISTICS 


PLANT OPERATIONAL 
CHARACTERISTICS 
& FINANCIAL FACTORS 



Figure B-3. Structure of Power Plant Simulation Computer Program 
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Capital costs as determined by the capital cost subprogram and 
capacity factor from the performance subprogram are inputs for the 
plant energy costs subprogram* The capacity factor is the ratio of 
the actual energy delivered throughout a year of operation to the 
energy that the plant could have delivered it if had operated 
continuously at rated plant power during the year* Other required 
inputs to the energy cost program include operational characteristics 
such as O&M costs and plant life. Financial factors, including 
considerations such as escalation rates, taxes, and insurance as 
presented in the previous section, are also needed as inputs to the 
energy cost subprogram* 

As seen from Figure B-3, the computer program contains an 
optimization feedback loop which determines minimum energy costs as a 
function of collector field area and storage capacity. A collector 
field area is first chosen and the storage size is increased until the 
minimum energy cost is reached. Since energy in excess of plant 
rating must be either wasted or stored according to the adopted 
operational ground rules, increasing storage reduces waste energy 
which, in turn, increases system efficiency and capacity factor. This 
reduces energy cost until the storage size reaches a value where 
further increases are no longer effective in reducing waste energy. 
After this point is reached, the added cost of further storage 
increases the energy cost. 

There is a particular capacity factor at the storage size 
corresponding to minimum energy cost. Selection of another collector 
area and a repetition of the process of increasing storage capacity, 
will provide a different minimum energy cost at a different capacity 
factor. By systematically varying collector field area and repeating 
the process, a curve of minimum energy costs as a function of capacity 
factor is determined. This constitutes the primary output of the 
computer power plant simulation program. 


B. MINIMUM ENERGY COST ENVELOPES 

Minimum energy cost envelopes are presented in Figures B~4, B-5, 
and B-6 for the paraboloidal dish, central receiver, and parabolic 
trough, respectively. For this set of three charts, envelope curves 
for five sites representative of major regions in the Continental 
United States are shown. The selected sites are: 




Region 

• 

Barstow, CA 

Southwest 

• 

Omaha , NE 

Midwest 

• 

Miami, FL 

Southeast 

• 

Maynard, MA 

Northeast 

• 

Medford, OR 

Northwest 
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Figure B-5. Minimum Energy Cost Envelopes for the Central 
Receiver System 
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Additionally, Albuquerque, NM, ia shown on Figure B-4 to 
illustrate that significant site-specific differences exist within any 
region. Both Albuquerque and Barstow are in the Southwest, but Barstow 
has a higher level of annual direct insolation and consequently 
achieves lower energy costs. 

The envelope curves are generally characterized by a sharp dip 
in energy costs at a low capacity factor and a rapid rise at high 
capacity factors. The dip at low capacity factors occurs when storage 
is no longer beneficial. The cost of a storage system is composed of 
a term proportional to plant power rating (or the rate at which energy 
is stored and extracted) and a term proportional to storage capacity 
(see Ref B-25). As long as even a small amount of energy is to be 
stored, the power-dependent cost term remains virtually unchanged for 
a fixed power from storage. Thus, at a low capacity factor when the 
storage system is no longer needed, the power-dependent capital cost 
is deleted and a corresponding decrease in energy costs results. 

As the average annual insolation decreases, the no-storage dip 
occurs at lower capacity factors, e.g., the dip occurs at a lower 
capacity factor at Maynard, MA, than at Barstow, CA. The annual total 
direct insolation at Maynard is 1253 kWh/m^ yr as compared to 2848 
kWh/m2 yr in Barstow, Regions of low insolation require a larger 
collector field area at a given capacity factor (or annual level of 
energy delivery). In these regions, there are some days where 
insolation intensity is comparable to that in high insolation 
regions. For these sunny days, the large field area results in excess 
energy that needs to be stored. Hence, storage is generally employed 
at lower capacity factors for low-insolation regions. 

The dip in energy cost is more pronounced when it occurs at low 
capacity factors. This results because the energy cost savings is 
proportional to the capital cost change (which is essentially 
constant) divided by the capacity factor (which is a measure of energy 
delivered). Therefore, when capacity factor is low, larger energy 
cost changes occur. 

As the collector area and capacity factor decrease beyond the 
dip, energy costs increase. This occurs because of the selected 
operational ground rules. The power systems were not allowed to 
operate at below ten percent of rated power. For plants of smaller 
collector area, this cut-off becomes increasingly significant. 

The minimum point in the dip corresponds approximately to the no- 
storage case where the plant is able to produce rated power only on 
those days having high solar intensities. For capacity factors lower 
than this minimum, the plant cannot produce rated power. This region 
is not of interest and is presented only in the context of identifying 
the minimum-energy-cost, no-storage plant. 

The rapid increase in the minimum energy cost envelope at high 
capacity factors is governed by characteristics of the insolation 
patterns. The maximum capacity factor that can be achieved is 
determined by the scheduled and unscheduled maintenance downtime. (A 
7 percent downtime is used for all plants in this study.) As this 
maximum level is approached, the plant must operate on a nearly 
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continuous basis. To operate on a continuous basis, large storage 
times are needed since, in all regions, there are periods of inclement 
weather where days of storage are required for the plant to produce 
power. To provide days of storage requires large collector field 
areas. The combination of large field areas and large storage results 
in high capital costs for a small increase in capacity factor and 
results in the rapid increase in energy costs at high capacity factors 
shown in Figures B-4, B-5, and B-6. 

At Medford, OR, the rapid increase in energy cost with 
increasing capacity factor occurs at a relatively low capacity 
factor. This effect is due to a large seasonal shift in insolation 
characteristics as shown in Appendix A. The area and storage 
combination that is best for one season performs poorly during the 
other seasons. This mismatch results in the need for large areas and 
storage capacities at unusually low capacity factors. 

All the minimum energy cost envelopes for the three systems 
exhibit the same basic trends. The lowest energy costs are achieved 
in the solar-intensive Southwest. At the lower insolation levels 
corresponding to other sites, the curves shift to higher energy 

costs. The two-axis tracking central receiver and paraboloidal dish 

achieve energy costs of about 50 mills/kWeh in the Southwest sunbelt, 

whereas the one-axis tracking parabolic trough has a value of about 75 

mills/kWe hr. The non-tracking CPC exhibits similar trends but at a 
higher energy costs. 


C. COLLECTOR FIELD AREAS AND STORAGE TIMES 

Each point along the minimum envelope curves in each of Figures 
B-4, B-5, and B-6, corresponds to a different variant of the power 
plant type presented by the figure. Each variant plant has a 
different combination of field area and storage time. If the capacity 
factor and regional location (insolation characteristics) are 
specified, the collector field area and storage time corresponding to 
a plant yielding minimum energy costs can be determined. Figures B-7 
through B-9 address collector field areas whereas Figures B-10 through 
B-12 pertain to storage. 

Collector areas for minimum energy cost plants are presented i- 
Figures B-7, B-8, and B-9 for capacity factors of 0.55, 0.40, and no 
storage (capacity factor variable), respectively. The following 
overall trends are shown by these charts: (1) Higher capacity factors 

require larger areas. (2) Regions of lower insolation need larger 
collector areas for any given capacity factor. (3) Two-axis tracking 
systems (PDS and CR) generally require smaller collector areas than 
one-axis tracking systems as represented by the parabolic trough 
(PT). The non-tracking CPC (now shown) required a larger collector 
area than the PT. 

At a capacity factor of 0.55 (Figure B-7), the areas required at 
Medford, OR are considerably higher than is indicated by the general 
trend. As previously described, this is due to the strong seasonal 
variation in insolation at Medford. At a capacity factor of 0.40 
(Figure »-S), ♦■he Medford data are consistent with general trends. 
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COLLECTOR AREA IN THOUSANDS OF SQUARE METERS 



BARSTOW, 
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Figure B-8. Effect of Regional Insolation Differences on 

Collector Field Area fur Capacity Factor = 0.40 
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Figure B-9. 


Effect of Regional Insolation Differences on 
Collector Field Area for No Storage 
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STORAGE TIME IN HOURS 
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Figure B-10. Effect of Regional Insolation D Terences on 
Storage Time for Capacity Factor - 0.55 


1 





Figure B— 11. Effect of Regional Insolation Differences on 
Storage Time for Capacity Factor ■ 0.4 
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CAPACITY FACTOR 







The no-storage case (Figure B-9) is presented primarily to show 
that relatively small areas are sufficient for all regions. These 
areas were estimated by interpolation in the dip region of the minimum 
energy cost enevelope curves, and this accounts for the scatter in 
some of these data. Detailed computer calculations to improve the 
accuracy of these estimates were not undertaken in this study since 
such improvement would not alter the study conclusions. 

The storage times corresponding to optimum plants at capacity 
factors of 0.55 and 0.40 (Figures B-7 and B-8) are presented in 
Figures B-10 and B-ll, respectively. General trends are that storage 
time increases with capacity factor and decreases with increasing 
insolation, i.e., the lowest storage requirements are in the sunbelt 
region. For the no-storage case, capacity factor is presented as a 
function of insolation level in Figure B-12. As discussed earlier, 
the capacity factor corresponding to the no-storage case decreases as 
the insolation level decreases. The points shown on Figure B-10 
through B-12 were determined by approximate interpolation procedures 
using the minimum-energy-cost envelopes. Although the procedure 
results in some data scatter, the general trends are well delineated. 
Interpolation accuracy could have been improved by using a fine mesh 
of collector areas and storage times to determine the envelope 
curves. However, this would have required additional costly computer 
runs, and improved accuracy was not considered to be warranted by the 
purposes of the present study. 


D. CAPITAL COST BREAKDOWN FOR SELECTED SYSTEMS 

Capital costs for the four selected systems are presented in 
Figure B-13 for three locations (Barstow, Omaha, and Maynard). The 
capacity factor is 0.55 for all the plants shown in this figure. The 
capital costs for the two-axis tracking systems (PDS and CR) are lower 
than the one-axis tracking PT and the non-tracking CPC at all 
locations. Lowest costs occur in the sunbelt region as represented by 
Barstow. 

As the insolation level decreases, progressively higher capital 
costs are required. This follows from the increased requirements for 
both collector field area and storage time as shown in Figures B-7 
through B-12. Capital costs at Maynard are roughly double those at 
Barstow. 

For each of the plants, capital costs are broken down into major 
cost categories. Collectors are the largest single cost item of the 
major categories shown. Electrical transport costs for PDS and 
thermal transport within the tower of CR are relatively small. For 
the CR, the transport costs, as well as the tower cost, are included 
in the collector cost. The pipeline network thermal transport systems 
used for both the PT and the CPC comprise a significant fraction of 
their total capital costs. 

Storage cost increases for locations such as Maynard relative to 
Barstow are directly related to the need for greater storage times, 
cf. Figure B-10. Also, les9 efficient systems, such as the CPC, 
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Figure B-13. Capital Costs for Selected Power Plants at 
Capacity Factor = 0.55 








require larger size thermal storage units since a relatively smaller 
fraction of the stored energy is converted to electricity. 

The energy conversion costs (including balance of plant) are the 
same for all the three systems which employ a central 10 MWe power 
generation unit (CR, PT, and CPC) since the same unit cost was assumed 
for all engine types. For the PDS system, conversion costs are linked 
to collector field area. The basic power generation nodule is a dish 
with a small heat engine/generator assembly at the focal point of each 
collector and, hence, as the collector field area is increased by 
adding more modules, the number of engines increases proportionately. 

Other capital includes indirect costs such as architect and 
engineering (A&E) and construction management fees. Spares and 
contingencies are also included in other capital. A value of fifteen 
percent of total direct costs is used for all systems. 
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I. INTRODUCTION 


The nominal performance and cost values used in this study 
anticipate the successful development of advanced technologies and a 
level of market penetration for solar systems that will result in 
high-volume mas a production of critical cost components such as 
concentrators. These projections were made in the context of the 
1990-2000 timeframe, and, accordingly, computer simulations (Appendix 
B) were based on plant startup in the year 2000. 

Variations in the extent to which advanced technologies are 
developed and the level of market penetration will result in 
variations about the selected nominal values. To delineate the effect 
of variation about the nominal projected values, selected parametric 
tradeoff and sensitivity analyses have been performed, and the 
corresponding results are presented in this appendix. 


II. BASIC TRADEOFFS 


For the different types of solar thermal power plants, the key 
factors governing comparative techno-economics are system efficiency 
and collector costs. This basic tradeoff and its ramifications are 
first described for four types of plants consisting of the (l) 
paraboloidal dish Stirling (PDS), (2) central receiver (CR), (3) 
parabolic trough (FT), and (4) Compound Parabolic Concentrator (CPC). 
Then, fundamental design tradeoffs are discussed in the context of the 
effect of site-specific or regional insolation differences. 


A. EFFICIENCY VS COST 

The effect of system design efficiency on busbar energy costs of 
the four selected advanced systems is presented in Figure C-l. System 
design efficiency is the ratio of delivered electrical energy to the 
design value chosen for the product of direct insolation and collector 
field area. For a fixed plant rating and specified annual energy 
delivery as denoted by capacity factor, power plant types having 
higher system design efficiencies require smaller collector field 
areas. 


Generally, more efficient plant types (CR and PDS) employ 
two-axis tracking collector configurations that are more costly per 
unit area than simpler one-axis tracking (PT) and non-tracking 
arrangements (CPC). Although lower collector costs were projected for 
the simpler collector systems (PT and CPC), the higher system 
efficiencies of the two-axis tracking systems (CR and PDS) exert a 
dominant influence which results in lower busbar costs for the 
systems. As illustrated by the charts for Barstow, CA, Omaha, NE, and 
Maynard, MA, this basic trend holds for all the regions examined in 
the study. 
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The effects of variations in concentrator costs and system 
design efficiency are delineated on Figure C-l. For the PDS system, 
the concentrator cost target range is 70-100 $/m?. A nominal value 
of 85 $/ffl2 was used in this study. The target range of concentrator 
cost represents a +18% vertical variation about the nominal value. 
This variation was applied to all the plant types to determine the 
energy cost range shown on Figure C-l. The cost values are given 
below: 


Concentrator Costs, $Aj2 


Type of Plant 

Low 

Nominal 

High 

PDS 

70 

8 r 

100 

CR 

53 

65 

77 

pt(D 

57 

70 

83 

cpcQ) 

41 

50 

59 


(l)l 


ncludes receiver costs. 


The nominal system efficiency was based on achieving 70% of 
Carnot. Present engines achieve about 60? of Carnot, whereas very 
advanced engine concepts are projected to attain values approaching 
80% of Carnot (Ref. C-l). A range of system design efficiencies 
corresponding to energy conversion efficiencies ranging from 60%, to 
80% of Carnot has been applied in Figure C-l and results in the 
horizontal variation about the nominal value. 


B. SITE-SPECIFIC INSOLATION EFFECTS ON POWER PLANT DESIGN 

For each of the power plant types considered, design tradeoffs 
can be performed to optimize plant design for a given set of 
site-specific factors. Insolation differences from site to site 
require detailed design perturbations in order to achieve optimum 
performance at each of the sites. Application of these extensive 
design optimization procedures was not within the scope of the present 
study. However, it is important to identify key tradeoffs and note, 
as a caveat, that some tradeoffs remain to be identified before 
further design improvements/refinements can be made to the power 
plants analyzed in this study. 

1 . Concentration Ratio and Operating Temperature 

The concentration ratio and operating temperature of the 
receiver can be tailored to match prevailing insolation 
characteristics with the objective of maximizing receiver efficiency. 
For a given operating temperature, receiver heat losses (reradiation, 
convection, and conduction) are essentially constant. Since receiver 
efficiency is the heat input rate minus the loss rate divided by the 
heat input rate, it follows that receiver efficiency improves as the 
heat input rate increases (for constant temperature). Thus, for a 
fixed collector geometry and given operating temperature, h ; ;hest 
receiver efficiencies are attained in the sunbelt regions having the 
highest insolation and, hence, the highest receiver heat input rate. 
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For sites where insolation is low, the reduction in receiver 
efficiency can be at least partially offset by allowing the receiver 
to operate at a lower temperature. However, reducing the temperature 
decreases the power conversion efficiency. For an optimum system, the 
product of these two efficiencies must be maximized. The systems 
treated in the present study were designed for optimal performance in 
the Southwest sunbelt. Thus, in other regions (particularly those 
having low insolation levels), a reduction in operating temperature 
might result in improved performance. 

Another variable that affects receiver performance is 
concentration ratio. Since receiver performance at a given 
temperature improves as the heat input rate increases, higher 
concentration ratios (which increase the heat input rate) are 
generally desirable. This is particularly true at high operating 
temperatures. Increases in concentration ratio require more accurate 
concentrator surfaces. Since the advanced projections are already 
based on the achievement of highly accurate reflective surfaces, a 
further improvement in accuracy to increase concentration ratio was 
not considered to be a feasible way of compensating for lower 
insolation levels. 

The CPC concept involves a further tradeoff. A part of the 
diffuse insolation is used to heat the receiver. As the concentration 
ratio is reduced, a greater portion of the diffuse insolation can 
reach the receiver because side-wall reflective surfaces are smaller 
and impose less limitations on the collector’s field of view. 

However, the direct beam insolation will be less concentrated. This 
results in a reduced heat rate to the receiver from beam insolation. 
Thus, concentration ratios for the CPC are chosen to achieve a balance 
between these effects so as to provide the highest overall performance. 

Argonne National Laboratory, developer of the CPC concept, 
recommends concentration ratios of 3 to 5 for all regions. For most 
regions, the diffuse component is very small; even for locations such 
as Miami, FL, which ts known to have a high level of diffuse 
insolation, the diffuse component reaches a value of only about 30 
percent of the total horizontal insolation. Thus, insolation is 
available primarily in the direct or beam form, and large departures 
from the concentration ratios deemed appropriate for sunbelt regions 
are expected to degrade the performance of solar thermal collectors. 

2. Engine Selection and Sizing 

Regional insolation levels can affect both engine selection and 
sizing. If the receiver ooerating temperature is decreased for 
regions of low insolation, the choice of engine type could he 
affected. For example, steam Rankine -ngines appear to be preferable 
to Brayton engines at 1000°F. Braytci* ngines provide a relatively 
higher performance at temperature above 1500°F where present steam 
Rankine engines and their advanced derivatives are unable to operate. 
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For systems employing thermal storage (CR, FT, and CPC), engine 
siting is not affected by insolation patterns. The engines are sized 
in accord with the desired plant rating. (A rating of 10 MWe was 
selected for this study). For the PDS system, which uses battery 
storage, the selection of engine size for the study was based on 
choosing engine rated power so as to accept peak levels of insolation 
available in the Southwest sunbelt regions. As seen from the 
histograms presented in Appendix A, most sites achieve peak insolation 
values approaching those of the sunbelt but, in many cases, for only a 
very small fraction of the total time during which insolation is 
received. Under these conditions (which occur primarily in regions of 
low insolation), the Stirling engine operates at part-load ratings 
(Appendix B) where efficiency is appreciably degraded. 

In these regions of low insolation, if a smaller engine of 
lesser power rating were employed, it would operate closer to its 
design efficiency during those periods when the bulk of the power was 
being generated. A higher system efficiency would result *rom this 
small-engine design, even though s small amount of power a: peak 
intensities would be wasted (via, e.g., collector defocusing). 
Further, the smaller engine would require a smaller oapitil 
investment. Thus, the possibility exists for improving the 
techno-economics of the PDS system in low-insolation regions. 


SECTION I T I, SENSITIVITY TO UNCERTAINTIES IN PROJECTIONS 

The sensitivities of solar thermal plant energy costs to 
variations in the cost and efficiency of all major subsystems were 
determined in Ref. C-2. Although updated cost and efficiency values 
are used in the present study, the trends shown in Ref. C-2 are 
generally still valid. Selected sensitivity results are presented 
herein only to make specific points. 


A. CONCENTRATOR UNIT COSTS 

Concentrator unit costs are the most critical cost driver. The 
low-cost values used in this study are based on high-volurae mass 
production (Refs. C-3 and C-4). Sensitivity of busbar energy cost to 
concentrator unit cost is shown on Figure C-2. Curves for the four 
selected systems are presented at Barstow, Omaha, and Maynard. 

It is seen that the variation in energy cost with concentrator 
unit cost becomes more pronounced in regions of low insolation such as 
Maynard. For a given capacity factor (0.4 is used on Figure C-2), 
collector field areas are larger in regions of low insolation. Hence, 
energy costs are relatively more sensitive to concentrator unit costs 
in regions of low insolation. 

further, the more efficient systems (CR and PDS) require smaller 
collector areas and are therefore less sensitive than the PT and CPC 
systems to collector unit costs. 









B. ENERGY CONVERSION COSTS 

For large 10 MW central energy conversion systems, a cost of 60 
$/kWe was used as a lower-bound, mass production target. At present, 
engines in this size category are not mass produced and cost about 
$300/kWe. Since this low-bound value is considered to be a 
particularly difficult target for 10 MWe engines, sensitivity curves 
of incremental energy cost as a function of energy conversion cost are 
presented in Figure C-3. This curve pertains to all three systems 
(CR, PT, and CPC) which employ 10 MWe central engines. 

Incremental energy costs read from the curves are t 9 be added to 
the energy cost values determined for an engine cost of 60 $/kWe, 
e.g., for an engine cost of 140 $/kWe at a capacity factor of 0.4, an 
incremental energy cost of 3-4 mills/kWe hr is to be added. 



Figure C-3. Sensitivity to Energy Conversion Cost 
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C. OPERATIONS AND MAINTENANCE COSTS 

The costs for operation and maintenance (O&M) of advanced solar 
thermal plants are subject to a high degree of uncertainty. A nominal 
value of 1% of total capital costs is employed in this study per EPRI 
(Ref. C-5). Sensitivity of busbar energy costs to O&M costs are 
presented in Figure C-4 for the four selected power systems when 
located at Barstow, California and operating at a capacity factor of 
0.55. 


Since the two-axis tracking systems (CR and PDS) have lower 
capital costs, they also have lower O&M costs. If the O&M for these 
systems increased from 1Z to 3% of capital, energy costs increase by 
approximately 15 mills/kWe hr. For the same variation in O&M, energy 
cost increases are even greater for the PT and CPC systems. 


D. STORAGE SYSTEM COSTS 

Storage system costs comprise a significant part of total plant 
capital costs at capacity factors in the neighborhood of ^0.5 (see 
Figure B-13 of Appendix B). Since there are large uncertainties 
regarding the costs of advanced storage systems, sensitivities to 
these costs are presented in Figures C-5 and C-6 for thermal and 
battery storage systems, respectively. 

In Figure C-5, the energy cost associated with the capital cost 
of thermal storage is presented as a function of the combined storage 
throughput and energy conversion system efficiency for a range of unit 
thermal storage cost values. For illustrative purposes, a capacity 
factor of 0.55 and storage time of 6 hours are assumed. Here, unit 
thermal storage cost is defined as the capital cost of the total 
storage subsystem divided by its maximum storage capacity in kWeh. 

In determining the size of a thermal storage system that is to 
be used in a solar thermal power plant, it is noted that, for a 
specified production of electrical energy from stored thermal energy, 
the size of the thermal storage system depends on the product of 
storage throughput and energy conversion system efficiency (Ref. 

C-6). As this product or combined efficiency increases, a greater 
portion of stored thermal energy is converted to electricity and, 
hence, storage size decreases. Also, it will not be generally 
possible to withdraw all of the stored energy, and this inability will 
require additional oversizing. For this study, a required oversizing 
of about 5% was assumed. 

The reduction in thermal storage size with increasing efficiency 
results in lower capital costs and associated energy costs as shown by 
the curves of Figure C-5. More efficient systems such as the CR thus 
gain an advantage over less efficient systems, particularly for high 
unit storage costs. This advantage is partially offset since 
low-efficiency systems generally employ lower temperature thermal 
storage having lower unit costs (see Table 3 of the body of this 
report) . 
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BUSBAR ENERGY COST, millsAWhr 








COMBINED STORAGE & ENERGY CONVERSION SYSTEM EFFICIENCY, % 


Figure C-5. Sensitivity to Thermal Storage Costs 



C-12 




The unit thermal storage cost range used in this study ranged 
from 21 to 26 $/kWeh, where unit costs were assumed to increase with 
temperature. There is considerable uncertainty in storage costs since 
advanced systems being projected in this study are in an early stage 
of development. Recent activities at NASA/LeRC have established 
targets or goals for some advanced systems in the range of 
approximately 5 to 15 $/kWeh. If these targets are achieved, energy 
costs of the order of 10 mills/kWeh lower than the nominal values 
presented in the present study will result. 

In addition to cost uncertainties, advanced batteries also face 
uncertainties regarding life. If battery life is less thap the 30 
year plant life, energy costs increase as shown in Figure 06. These 
energy cost increases for battery replacements are determined via the 
present value economic methodology of Ref. 07. 

A unit battery storage cost range of 30 to 60 $/kWeh is used in 
Figure 06 to illustrate how unit costs affect the impact of battery 
replacements. For the selected unit cost range, energy costs increase 
from 5 to 10 mills/kWeh if battery life decreases from 30 to 15 years 
(corresponding to one replacement). At short battery lifetimes of the 
order of 5 years (representative of present lead-acid batteries) 
substantial cost increases occur. 

Thus, there is considerable incentive for developing an advanced 
battery, such as the Redox, which has potential for both low cost and 
long life. Other storage systems, such as underground pumped hydro 
and compressed air, can be used to store electro-mechanical energy, 
but, if cost and performance targets are achieved, battery systems are 
desirable because of implementation flexibility. 
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